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ABSTRACT 



The nature of galajcies selected at submillimeter wavelengths (SMGs, 5*850 ^ 3 mjy), some of the bolometrically most 
luminous objects at high redshifts, is still elusive. In particular their star formation histories and source of emission are 
not accurately constrained. In this paper we introduce a new approach to analyse the SMC data. Namely, we present 
the first self-consistent UV-to-radio spectral energy distribution fits of 76 SMGs with spectroscopic redshifts using all 
photometric datapoints from ultraviolet to radio simultaneously. We find that they are highly star-forming (median 
star formation rate 713 M© yr"'^ for SMGs at z > 0.5), moderately dust-obscured (median Av ^ 2 mag), hosting 
significant stellar populations (median stellar mass 3.7 X 10" Mq) of which only a minor part has been formed in the 
ongoing starburst episode. This implies that in the past, SMGs experienced either another starburst episode or merger 
with several galaxies. The properties of SMGs suggest that they are progenitors of present-day elliptical galaxies. We 
find that these bright SMGs contribute significantly to the cosmic star formation rate density (~ 20%) and stellar 
mass density (~ 30-50%) at redshifts 2-4. Using number counts at low fiuxes we find that as much as 80% of the 
cosmic star formation at these redshifts took place in SMGs brighter than 0.1 mJy. We find evidence that a linear 
infrared-radio correlation holds for SMGs in an unchanged form up to redshift of 3.6, though its normalization is offset 
from the local relation by a factor of ~ 2.1 towards higher radio luminosities. We present a compilation of photometry 
data of SMGs and determinations of cosmic SFR and stellar mass densities. 
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1. Introduction 



Submillimeter galaxies (SMGs; 
were discovered at 850 /im (5*850 
Submillimetre Common-User Bolometer Array (SCUBA 
[Holland et al.l . [T999i ) mounted on the James Clerk Maxwell 
Telescope (JCMT). Due to the coarse resolution of 
SCUBA, localizations derived from high-resolution radio 
maps had to be used t o mea sure their spectroscopic red- 
shifts (" Chapman et al.l [2005). Lots of studies have ad- 
dressed the issue of characteriz ing the nat ure of SMG s 
(lEgami et all. 120041: iGreve et al.l. [2004. 20 05t ISmail et al 

11 



galaxies: evolution - galaxies: high-redshift - galaxies: ISM - galaxies: starburst - sub- 



to solve several issues, including the characterization of the 
star formation histories of SMGs and their dominant source 
of emission. 



see iBlain et al.l I2002D 
> 3 mJy) by the 



20041; 'Swinb ank et al.L 12004 l2006l I200& Takagi et al J."l2004 ; 



Alexander et al .L 120051: iBorvs et al.l. 12005: Kovacs et al 



2006; Laurent et al., 2006; Pooe et al., 200 0: iTacconi et al 



2006. 



2008, 



2008i: ,Takata et al.,. ,2006i: ,Younger et al.l 



Coppin et al.l. 



.200931 IClements et al.l. l2008t 

Dve et al.', '2008', '2009'; 'Hainline', '2008'; 'Hainline et alj 
Pcrcra et al., 2008; Scott ct al., 2008 



2007 



2008; 



20091; 

Austcrnian n et al 
Devlin et all. 120091: lEales et al.l. 12 009: Murphv et al. 

L|, .2009: Weifi et al. 



2009; . ■ ■_ 

2009; Murphvl. 1200^ : 
2009b; Aravena et al 



iTamura et all . ,2009 
I2OIOI some of these works were 
based on surveys with sensitivity worse than 3 mJy quoted 
above). However they were usually based on limited sam- 
ples (< 20 sources), limited wavelength coverage or pho- 
tometric redshifts. These limitations have made it difhcult 



An important open question concerns the contribution 
of SMGs to cosmic stellar mass assembly. This is impor- 
tant, because in order to understand galaxy evolution, the 
build-up of stellar mass must be mapped out to high red- 
shifts. It is usually parametrized by the total star formation 
rat e (SFR) density per unit comoving volum e, (psfr; see 
e.g. lHopkinsll2004j:lHopkins fc Beacoinl . l2006D . At high red- 
shifts it is difficult to disentangle the contribution to psFR 
from galaxy populations of different masses due to incom- 
pleteness at low luminosities. 

Another approach to study stellar mass assembly is 
to consider directly the stellar mass density per unit co- 
moving volume, which is equivalent to the integrated 
PSFR over the age of the Universe. It is established that 

grows with cosmic time (stellar mass is accumulat- 
ing; Drorv ct al., 2005; Fontana ct al., 2006; Eisner et all 
[2OO8I: iPerez-Gonzalez et al.l . 120081: iMarchesini et al.l . l2009h . 
but the contribution from different galaxy popula- 
tions i s not well-deterniined. Spitzer observations of 
SMGs (Egami ct al.','2004VFravcr ct al.','2004;'lvison et aL, 
2004; Borvs ct al., 2005; Ashbv ct al ,20Q 6 ; Laurent et a", 
20 061: IPope et all. 120061: lOve et al.l . WM iHahdind . I2OOS ; 
iHainline et al.Ll2009f ) have enabled studies of the rest-frame 
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near-infrared (near-IR) part of the spectrum, where old 
stellar populations are dominant - an important step for- 
ward in getting full spectral energy distributions and ac- 
curate estimates of stellar masses of SMGs. The results in- 
dicate that SMGs are among the most massive galaxies in 
the Universe. 

The dominant source of emission from SMGs is dust re- 
processed emission either from young stars or active galac- 
tic nuclei (AGNs). One way to test it is to compare the 
infrared (IR) and radio luminosities of SMGs, because, at 
least locally, star-forming galaxies follow a rem arkably tight 
corre l ation bet ween IR and radio luminosities (jHelou et al.l , 
119851: iGondonL |I992). The correlation is believed to re- 
sult from the fact that both IR and radio emissions are 
related to short-lived massive stars: the former originates 
from dust heated by ultraviolet (UV) light from blue, mas- 
sive stars and the latter from synchrotron emission of elec- 
trons produced in supernova remnants. Therefore, a re- 
lation consistent with the local one is an indication of 
star formation dominating both the IR and radio emis- 



at redshifts z < 


1 (iGarrettl. 20021: iGruDoioni et al.l. 2003 


Aooleton et al. 




2004; iBovle et al., 120071: iMarleau et al. 


20071: Vlahakis et al.l. 120071 lYane et al.l. 12007). At higher 



redshifts sample sizes are small making it difficult t o draw 
robust conclusions ([Applcton et al.^ 2004; Ko vacs et al 



2006; B eswick et al.l. 2008: Ib ar et al I 2008; S aiina et al 
'2008i,Garn et al.l.l2009c.MurDhv et al.ri20q9:,MurDhvl.l2009 : 



Rieke et al.l . 120091: ISevmour et al.l . 120091: lYounger et al . . 
2009b'). Th e only sign of evolution was reported by 
Ivispn et al.l (|2009t ) based on stacking analysis of the 24 /im- 
selected galaxies, though possibly interpreted as a selection 
effect. 

The objective of this paper is to model for the first time 
the entire UV-to-radio spectral energy distributions of a 
statistically significant sample of SMGs in a self-consistent 
way. Using these models we i) consistently derive the prop- 
erties of SMGs using all available data to characterize their 
nature and determine the dominant emission mechanism; 
it) estimate the contribution of SMGs to the cosmic SFR 
and stellar mass densities; iii) investigate whether the lo- 
cal IR-radio correlation holds at high redshifts in an un- 
changed form. In Section [2] our SMG sample is presented. 
Our methodology is outlined in Section [31 We derive the 
properties of SMGs in Section 2) and discuss the implica- 
tions in Section [5] Section |6] closes with our conclusions. 
We use a cosmological model with Hq = 70 km s^^ Mpc~^, 

= 0.7 and fi,™ = 0.3. 

2. Sample 

We base our analysis on 76 SMG s {Ss5o ^ 3 mJy) from the 
sample of lChapman et all (|2005[ ). all with spectroscopically 
measured redshifts spanning a range of 0.080-3.623. 

The way the sample is selected involves complex bi- 
ases, which are difficult to fully quantify and account for. 
The parent sample of Chapman et al. (2005) consists of 150 
SMGs out of which 104 have radio identifications. The sam- 
ple discussed here (76 galaxies) consists of the SMGs for 
which redshifts have been measured (spectroscopic com- 
pleteness ^ 75%). All this implies that the sample is biased 
against: i) faint submillimeter emitters (low dust content 
and/or hot dust, influence mostly the low-z portion of the 
sample); ii) faint radio emitters (high-z and cold dust, see 



Figure 3 of IChapman et al.l . l2005|) ; iii) faint optical emit- 
ters (difficult to obtain spectra); iv) z ~ 1.2-1.8 ("redshift 
desert" where no emission lines enter the observable wave- 
lengths). At low redshifts {z < 1) the sample may also 
be incomplete due to a limited sky area (and therefore - 
volume) coverage making it difficult to detect rare strong 
submillimeter emitters (for det ails on the SMG selection ef- 
fects see also Fi gure 2 of Blain et al.l (f2b04i) and discussion 



in Section 4.4 of iMichalowski et al.l (|2008D ) 



It is important to estimate what the infiuence of these 
selection effects on our results is. In total we analyse ~ 
50% (76/150) of the parent sample. Additionally, 25 radio- 
detected SMGs without spectroscopic redshifts have sim- 
ilar long-wavelength pro perties compa r ed to the redshift 
sample (see Figure 1 of iChapman et al] . l2005 l). so their ab- 
sence from the sample probably does not significantly bias 
our results. The same is true for the SMGs in the "red- 
shift desert" , since they are missed not due to their inher- 
ent properties. The remaining 46 radio-nondetected SMGs 
(~ 30%) could in principle have very different properties 
than our sample resulting in a potential limitation in our 
analysis. 

Even if most of the SMGs without spectroscopic red- 
shifts are similar to those in our sample, the incomplete- 
ness at z < 1.8 implies that the estimates of SMG densities 
(Sections 15.3.11 15.3.21 and 15.2.31) in the three low-redshift 
bins (see Sectioning) are strict lower limits. 

Due to the negative K-correction at submillimeter wave- 
lengths, SMGs at z > 0.5 form a sample with homoge - 
nous IR luminosity (jBlain fc Longaiii 119961: iBlainl . [1997). 
However, SCUBA sources at z < 0.5 belong to a different 
population of objects and are intrinsically fainter. The lim- 
ited volume coverage at these low redshifts makes the sam- 
ple of these objects small and incomplete. This prevents 
a separate study of their properties. We did not take into 
account these sources when we computed median values of 
the properties of SMGs. 

The photometric datapoints (Tables lA.ll a nd | A.2lin ap- 
pendisQ) were collect ed from the literature: llvison et al.l 
(20 02i IK, ra dio), [ivison e t al\ (f2005l R. 1. 2 mm) 

mo. 



Chapman et al. I (l2003k Vn. IChapman "^nn (l2005l 



BR, 850 /xm, r adioV [ Capak et al.r(l2004 UBVB Iz HK) 



IClements et all l2004 U BVIK),Es rami et al. 



24 iim) iGreve etal 
UK) 



T2OO4 1 .2 mm), iSmail et al 

iFomalont et al.l (I2006L R z) 



1 



(2004, 



2004 . 



2006 



350 Mm, 1.2 m m 
iTacconi et al 



iKovacs et ah. ., 

.Laurent et al ] (I2006r 350 /xm , 1.1 m m) , 
20061 1 3 ^i^. lPope etal.1 (l2006l R 



24 ^m), iHuvnh et all (|2007l 160 um). iHainlind (|2008l 
3.6,4.8,5.6,8.0,24, 70 /xm). We have n ot used the exist- 
ing mid-IR spectra (IValiante et al.. 20071: iPope et al.l . 120081: 



iMenendez-Delmestre et aL . 12007. ,2009i ) , but for complete- 



ness we have indicated in Table ET] those SMGs for which 
Spitzer/TRS spectra exist. 

3. Methodology 

3.1. SED modeling 

In order to model the spectral energy distributions (SEDs) 
of SMGs, we use all the photometric datapoints simultane- 
ously. This has the advantage that all the galaxy properties 



^ For convenience we make the compilation available in elec- 
tronic form. We suggest that the original data source be con- 
sulted and referred to appropriately. 




10"^ 10"^ 10° 10^ 10^ 10^ 10"^ 10^ 10^ 
Rest Frame Wavelength (|Lim) 



Fig. 1 Median spectral energy distribution (SED) of SMGs {thick lines) and SEDs of individual SMGs {thin lines). Dotted 
lines indicate z < 0.5 objects. Shaded areas enclose 90% of the SEDs. Top: all SEDs were divided by the corresponding 
850 datapoint and scaled, so that the median SED has a flux of 5 mJy at the rest-frame 283 /zm (observed 850 /Ltm at 
z ~ 2). Bottom: SEDs were normalized to an infrared star formation rate of 100 Mq yr"-"^. 



4 



Michalowski et al.: Evolution of submillimeter galaxies 




Fig. 2 Redshift evolutio n of the proper t ies (f ull circles, see Table IA.3I in appendix) of the sample of 76 SMGs with 
spectroscopic redshifts (jChapman et all I2005D . Small symbols indicate z < 0.5 objects. Typical errors (Section H]) are 
shown as crosses. From top-left to bottom-right: star formation rate (SFR) derived from spectral energy distribution 
modeling, ultraviolet, infrared and radio emission, SFR per unit stellar mass (= SFRir/A/»), stellar mass, fraction of 
stellar population formed during the ongoing starburst, stellar mass-to-light ratio, dust mass and temperature. In the 
SFRiR panel, we also show the minimum average SFRs (see Section I5.2.ip required to build up the total stellar mass 
within the age of the Universe at a given redshift {empty circles) and to build up the fraction of stellar population that 
was not formed during the ongoing starburst {plus signs). The location of plus signs indicates that SMGs must have 
been highly star-forming even before the onset of the ongoing starburst. When empty circles and plus signs overlap, the 
contribution of the ongoing starburst to the total stellar mass of a galaxy is negligible (i.e. Mburst/-^* ~ 0). 
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are derived consistently regardless of the wavelength regime 
in which those properties shape the SEDs (for example, re- 
cent star formation governs the UV and far-IR parts of a 
spectrum of a galaxy, whereas accumulated stellar mass is 
responsible for near-IR emission) . Moreover in the full SED 
modeling no single datapoint drives the fit alone. 

We utilized the set of 35 000 models from 
llglesias-Paramo et alj (I2007D developed in GRASIL 
( Silva et al ■L 1199811^1 based on numerical calculations of 
radiative transfer within a galaxy. They cover a broad 
range of galaxy properties from quiescent to starburst. 
Their star formation histories are assumed to be a smooth 
Schmidt-type law (SFR proportiona l to the gas mass 
to some power, see ISilva et al.l . Il998l for details) with a 
starburst (if any) on top of that starting 50 Myr before 
the time of the evolution of a galaxy at which the SED 
is computed. Addit ionally we fitt ed templates based on 
nearby galaxies ( Silva et al.l. 1998| ) and gamma-ray burst 
host galaxies (1 Michalowski et aTT 120081 ) . We simultane- 
ously used all the photometric datapoints from UV to 
radio (Tables lA.ll and IA.2P . In cases where the data given 
by different authors were contradictory, we disregarded the 
obvious outliers. We scaled the SEDs to match the data 
and chose the one with the lowest x^. 

Based on the best fits we derived t he pro perties of the 
galaxies as explained in iMicha low ski et al.l |2008, 2009). 
In particular, SFRs, stellar (M*) and starburst (Mburst) 
masses were given as output from GRASIL, rest- frame UV 
and K {Lk) monochromatic luminosities were interpolated 
from the best-fitting SEDs, whereas IR luminosities (Lir) 
were integrated in a range 8-1000 /im, UV and IR SFRs 
(SFRiR was adopted for all subsequent calculations, be- 
cause SFRuv is on avera ge two orders of magnitude lower) 
were calculated using KennicuttI (|1998l ). dust masses (M^) 
were calcu lated from the 85 Mm d etections using equa- 
tion (5) of IMichalowski et al.l (|2009D and radio SFRs were 
calculated from the 20 cm detections using the empirical 
formula of'Bel (|2003D (see Section 4.2 of lMichalowski et al.l . 
[2009). Dust temperatures (T^) were estimated by identify- 
ing the peak of the dust emission and assuming an emissiv- 
ity index /3 = 1.3. The average extinction in the rest-frame 
y-band was calculated from the unextinguished starlight 
given in GRASIL: Ay = 2.5 log(unextinguished y-band 
starlight / observed y-band starlight). IR-radio correla- 
tion parameters were calculated according to the formula 
q = log(LiR[L0]/3.75 X lOiVL,i.4GHz[i0Hz-i]), where 
i'j/i.4GHz is a rest-frame 1.4 GHz luminosity density com- 
puted from the observed 1.4 GHz flux assuming a spectral 
slope of —0.75. 

3.2. Volume densities 

In order to calculate the SFR density, the stellar den- 
sity and the dust mass densities per unit comoving vol- 
ume, psFR, P* and pdust, we used the following an- 
gular areas for the submillimeter surveys (Table 1 of 
iGhanman etHI . ,200 5): CFRS-03: 60 arcmin^ and CFRS- 
14: 48 arcmin^ (jwebb et al.l . l2003b[ ). Lockman Hole: 12 2 
arcmin^ and ELAIS- N2: 102 arcmin^ (IScott et all l2002l) . 
HDF-N: 100 arcmin^ (IChaoman et al.L 2001). SSA-13 and 
SSA-22: 100 arcmin^ each (jGhapman et al.l l2003ah . total- 
ing 632 arcmin^. 



We divided our sample into four high-redshift bins 
(Table [1]) with approximately the same number of SMGs 
plus an additional bin for z < 0.5 sources (see Section 
12]). The densities in each bin were calculated as a sum 
of SFRiR (or M*, or M^) of ah SMGs in this bin di- 
vided by its comoving volume (a similar approach to cal- 
culate the SF R and number volume densities o f SMGs 
was taken by ICoppin et al.l. 120091: iDaddi et al.l l2009bl: 
lYounger et al.l l2009al: IWang et al.l . l2009D . The volume den- 
sities (Column 2) were found using the total area from the 
previous paragraph. 

We removed the contribution of ten SMG^, which were 
observed by SCUBA in the photometry mode (as opposed 
to the blank-fie ld mapping mode) targ eting optically-faint 
radio galaxies (jGhapman et al.l . 12005 ). These objects fall 
outside the fields discussed here. 

The method is therefore to analyse the fraction of the 
sky observed by SCUBA and estimate the number of SMGs 
and their volume densities. However, the true number of 
SMGs in our fields could be higher. On the other hand, re- 
gardless of the selection effects, the true number of SMGs in 
our fields cannot be lower than the number of SMGs in our 
sample. In turn, the true values of SFR and densities 
cannot be lower than the values we derive. Therefore our 
results on volume densities should be regarded as robust 
lower limits. 

H aving this in mind we note that the parent sam- 
ple of ICha pman et al.l (|2005D includes only 29% of all the 
SMGs detected in the used survey fields (compare with 
IScott et all . 120021: IWebb et al.l . l2003bl lah . Therefore even if 
we analysed the full parent sample the estimated densities 
would be conservative lower limits. We attempt to correct 
for this incompleteness by assuming that the parent sam- 
ple of Chapman et al. (2005,) is a fair representation of the 
total population. In this case our numbers should be multi- 
plied by 3.5 (~ 1/29%). This correction should in principle 
be derived separately for each redshift bin, but the missing 
redshift information for the majority of the SMGs in the 
used survey fields makes such calculation impossible. We 
note that this correction does not remove the bias against 
SMGs that are faint at radio and optical wavelengths, as 
discussed in Section [2l 

We have not applied a volume density correction for 
the AGN contribution, because it is at most minor. Even 
though a fraction of SMGs host AGNs and a few in- 
dividual SMGs have been shown to exhibit a signifi- 
cant AGN contribution to their emission, it is established 
that on average AGN activity is responsible for at most 
^ 10-20% of th e bolometric infrared emission of SMG s 
(lAlexander et all. l2005l 120081: iMenendez-Delmestre et al 
2007', '2009^ Valiante et al.', l2007t [Pope et al.', '2008|: 
Hainlinc et aL 



7t [Pop 

2009; Murphy et all [2009; Watabc et al 
'20Q9D. Therefore a potential error associated with the AGN 
contribution in our analysis of a statistically significant 
sample is smaller than the systemati c uncertainty (e.g . 30% 
error of luminosity-SFR conversion; iKenniciItfl . 



The percentage contribution of SMGs to the SFR and 
M* densities (Columns 5 and 7 of Table[T]) was calculated as 



|http: //adlibitum. oat . ts . astro . it/silva/def ault .htmll 



^ SMMJ123553.26-h621337.7, 
SMMJ123600.10-f 620253. 5, 
SMMJ123606.85-f 621021.4, 
SMMJ163706.51+405313.8, 
SMMJ221806.77+001245.7 



SMMJ123555.14-f 620901. 7, 
SMMJ123600. 15+621047.2, 
SMMJ123716.01+620323.3, 
SMMJ221804.42-H002154.4, 
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PSMo/ifSMG + Pother), where psMG is the density of SMGs 
at each redshift bin (Columns 4 and 6) and Pothor is the den- 
sity of other galaxies assumed to be an average of determi- 
nations (excluding lower limits) reported by other authors 
(Figure |31 Tables IA.4I and IA.5I in appendix) , for which the 
redshift ranges overlap with our bins. This way of calculat- 
ing the contribution is justified if SMGs do not enter the 
"other" samples of galaxies. This is usually the case because 
SMGs are faint in the optical. However, if this was not ful- 
filled, the real percentage contribution of SMGs would be 
even higher. 



4. Results 

The best fitf0 are shown in Figure ETD and the median SEDs 
(in flux and luminosity domains) are shown in Figure [1] 

The resulting properties of the galaxies are listed in 
Table [ST3] and shown in Figure [2] as a function of redshift. 
We notice similar trends to Hainline, (2008) that lower- z 
SMGs are less luminous and colder (see her Figures 4.7 
and 4.9). 

In two cases we obtained much better fits using 



the templates of Silva et al.l (jl998[ ) instead of those of 
llglesias-Paramo et al.l (|2007[ ). namely, an HR 10 template 
for SMMJ105151.69-h572636.0 and a spiral Sc template 
for SMMJ221733.12-f001120.2. In 9 case^ where our fits 
strongly underpredict th e 850 fim datapo i nt we adopted the 
LiR and estimates of IChaoman et al] (|2005[ ). 

The determination of the IR luminosity suffers from sys- 
tematic uncertainties depending on the choice of the SED 
template. Our approach of using all the optical, submil- 
limeter and radio data to constrain the shape of the SED 
results in a moderate syst ematic error in t he IR luminosity 
(less tha n a fac tor of 2: lBen et all . 120071 ). The choice of a 
ISalpeteij (figsl IMF with cutoffs of 0.15 and 120 Mq intro- 
duces a maximum systematic error of a factor of ^ 2 in the 
determination of the stellar masses and SFRs (jErb et all . 
|2006[ ). iBell et al.l ()2007D have also found that random er- 
rors in stellar mass are less than a factor of 2. Estimates 
of dust temperatures have uncertainties of ~ 5-10 K domi- 
nated by the unknown value of the emissivity index, /3. The 
SFR determination based on radio observations is accurate 
up to 30% since it agrees with the detailed spec trophoto- 
metric SED fitting (jMichalowski fc Hiorthl . |2007[ ). The un- 
certainties in q (defi ned in Section I5TT|) are ~ 0.3 (see also 
iKovacs et al.l . l2006l) . dominated by the error in Ljr. 

In order to assess the infiuence of the choice of emissivity 
index /3 = 1.3 on the dust mass estimates, we recalculated 
the dust temperatures and masses in a range of (3 of 1-2. 
The resulting error was less than a factor of 3.5. 

This is illustrated on Figure [31 where we present a more 
systematic analysis of this problem. We calculated the dust 
mass of a mock galaxy with = 40 K (this choice does 
not infiuence the results) using (3 in the range 1-2 assum- 
ing a fiux density of 5 mJy at a variety of infrared rest- 
wavelengths probed by observations. Then we normalized 



The SED fits can be downloaded from 



Ihttp : //archive . dark- cosmology . dk 



" SMMJ030226. 17+000624.5, 
SMMJ030236.15-h000817.1, 
SMMJ123636. 75-^621156.1, 
SMMJ123721.87-h621035.3, 
SMMJ221724.69-I-001242.1 



SMMJ030231.81+001031.3, 
SMMJ030238.62-f 001106.3, 
SMMJ123651. 764-621221. 3, 
SMMJ163639.0H-405635.9, 




1.4 1.6 
Emissivity index p 



Fig. 3 Derived dust mass of a mock galaxy with dust tem- 
perature = 40 K and a fiux density of 5 mJy at sev- 
eral infrared rest-wavelengths as a function of the assumed 
emissivity index /3. For each wavelength the dust masses 
were normalized to 1 at /3 = 1.5. The spread of the derived 
dust masses shows that the uncertainty of the dust mass 
resulting from unknown /3 is a factor of a few. 

dust masses to 1 at /3 = 1.5. We conclude that as long as 
the observations probe wavelengths longer than ~ 150 fxni 
{z < 4.7 for observed wavelength of 850 /xm), then the error 
on the dust mass resulting from unknown (3 is less than a 
factor of 5. 

None of these errors significantly affects our conclusions, 
because the inferred nature of SMGs would not be differ- 
ent even in the worst case scenario when all systematic 
errors work in one direction (increasing or decreasing the 
obtained values). Moreover, we analyse a statistically sig- 
nificant sample of 76 galaxies, so random errors of a factor 
of 2 are reduced to < 20% when an error of a mean is 
considered. 

Table |T] contains the volume densities and mean IR- 
radio correlation parameter divided into five redshift bins 
(see Section 122]). The uncertainties quoted on psFR and 
include the systematic 30% uncertainty of the Lir to SFR 
conversion (jKennicuttl . Il998l ) an d a factor of ^ 2 sys t emati c 
uncertainty in the stellar mass (jMichalowski et al.l I2008D . 
The systematic error resulting from our incompleteness cor- 
rection (Section 13. 2p is likely a factor of a few. 

5. Discussion 

5.1. Spectral energy distributions of SMGs 

We have presented the first successful attempt to fit the 
entire UV-to-radio SEDs of SMGs in a self-consistent way 
taking into account all the available data simultaneously. 
Our study provides evidence that GRASIL models can re- 
produce the SMG data. Namely, we found good fits for all 
SMGs in our sample with the best IR/submillimeter wave- 
length coverage^ except of SMMJ105238.30-f 572435.8. 

As is evident from Figure [U regardless of whether SEDs 
were normalized to the same observed 850 /im datapoint 



" SMMJ105201. 25-1-572445. 7, 
SMMJ163650.43+405734.5, 
SMMJ163706.51+405313.8 



SMMJ105230.73+572209.5, 
SMMJ163658. 194-410523.8, 
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Table 1 Mean values for SMGs in redshift bins 





Volume 


log pm 


PSFR 




log p* 






log Pdust 


2 
(1) 


(10 Mpc ) 
(2) 


(Lq Mpc"^) 
(3) 


(Meyr Mpc ^) 
(4) 


7o 

(5) 


(Mq Mpc ) 
(6) 


% 

(7) 


q 

(8) 


(Mq Mpc ) 
(9) 


0.080 - 0.500 


0.12 


7Qn + U.li 
' •'^■^-0.14 


0.0018 ± 0.0005 


5ti 


"■'^'J-o.ie 


1+0.4 

-■--0.2 


2.54 ±0.12 




0.510 - 1.316 


1.03 


7 oi +0.06 

"•-'-^-o.os 


0.0111 ±0.0016 




'•^0-0.08 




2.52 ±0.06 


4 O2+0.04 


1.408 - 2.142 


1.63 


0.0228 ± 0.0027 


u+_l 


7 1 O+0.16 


11« 


2.29 ±0.08 


o 04+0.04 


2.148 - 2.565 


0.89 


8.45lH^ 


0.0486 ± 0.0054 




7 61+°-^2 

' -"^-0.08 


511?? 


2.27 ±0.09 


4 nc-+0.04 


2.578 - 3.623 


2.17 


8.30iH^ 


0.0341 ± 0.0040 


20+1 


7 r,o+0.12 
' -^"-O-OT 


31+^* 


2.25 ±0.08 


oe+0.03 
•3-OD-0.03 



Note. — Column (1): redshift range of the bins. Column (2): comoving volume of these bins (calculated in Section 13. 2p . 
Column (3): IR luminosty density of SMGs. Column (4): Resulting IR SFR density of SMGs (Section [53TT|) . Column (5): contri- 
bution of SMGs to the cosmic SFR density (calculated in Section [3. 2[l . Column (6): stellar mass density of SMGs (Section 15. 3. 2|l . 
Column (7): contribution of SMGs to the cosmic M, density (calculated in Section [3. 2 1) . Column (8): mean (and error of the mean) 
FlR-radio correlation parameter for SMGs (Section l5.4.1|l . Column (9); dust mass density of SMGs (Section 15.2.311 . Columns 3-7 
and 9 have been corrected for incompleteness by a factor of 3.5 (Section 13. 2p . 



or SFRiR, the scatter at optical and near-IR wavelengths 
is significant, showing that SMGs exhibit a wide range of 
stellar population properties (as also noted bv llvison et al.l . 
[2002). This implies the need for an SED template library 
in SMG studies, as opposed to single-template fitting. 

Having constrained the SEDs of SMGs we now turn to a 
discussion of what we can learn about these galaxies using 
the best-fitting models. 

5.2. Properties of SMGs 
5.2.1. Star formation rates 

The very high (current) SFRs of SMGs (median 713 M© 
yr~^, Column 5 of Table IA.3I and Figure ^ place them 
among the most powerful starburst galaxies in the Universe. 
Such extreme SFRs li kely result from r najor mergers (e.g . 



Chapman et al.', '2004; 'Swinbank et al l 120041: |G reve et all 



2005 ; Tacco ni et al., 2006, 2008; Youn ger et al.L 12007 2008: 
Berciano Alba et al.l . 120101: iNaravanan et al.l . l2009l 12010 ) 
and cannot be sustained for a long period (after a few 
hui idredMyr at most the gas reservoi r shou ld be depleted; 
see lGreve etall l2005t fHainline et aLl . l2006[ ). 

On the other hand, their extinction-uncorrected UV 
SFRs are two orders of magnitude lower (median 7 Mq 
yr"-'^, Column 4). This implies that the majority of star 
formation in SMGs is hidden by dust. Therefore, optical 
observations alone are not sufficient to investigate their na- 
ture and contribution to cosmic star formation. 

Using stellar masses of SMGs we placed lower lim- 
its on the time-averaged SFRs required to build their 
stellar masses within the age of the Universe (= 
M*/age of the Universe at given redshift), shown as empty 
circles on Figure Their median value of ~ 130 Mq yr~^ 
indicates that SMGs had to be relatively highly star- 
forming throughout the age of the Universe to build up 
their stellar populations at a constant rate. Even if our 
estimates of stellar masses were underestimated by a fac- 
tor of a few due to systematic uncertainties (Section H]), 
the SMGs would have had to be luminous infrared galaxies 
(LIRGs with SFR > 20 Mq yr~i) during their evolution. 

Having constrained the mass of stars formed during the 
ongoing starburst episode, Mbmst, we can further constrain 
the minimum average SFR of SMGs before the onset of this 



starburst, = (M* — Afburst)/age of the Universe (plus signs 
on Figured]). The median is still high, ^ 100 Mq yr"-'^, so 
SMGs must have been highly star-forming in the past too. 
At redshifts 2-3 the age of the Universe is ~ 3-2 Gyr and it 
is unlikely that a galaxy can sustain this high SFR over such 
a long period. Therefore we conclude that either the stellar 
masses of SMGs have been formed in at least two strong 
(> 100 Afg yr~^) starburst episodes or continuously over 
the period of 2-3 Gyr but in several smaller galaxies that 
eventually merged. In order to build up the stellar mass of 
one SMG, five such galaxies would need to form stars con- 
tinuously at a rate of 20 Mq yr~^, a value more likely to be 
sustainable over sever al Gyr. The latter scenario is consis- 
tent with the results of lDve et al.l (pOOS) based on observed 



optical to mid-IR data of 51 SMGs with photometric red- 
shifts. They found that approximately half the stellar mass 
in SMGs has been formed over a long (~ 1-2 Gyr) period of 
approximately constant star formation activity. The pos- 
sibility that a significant part of stellar mass in SMGs was 
for med before the o ngoing starburst has also been suggested 
bv iHainlind ([2003), who compared the build-up timescale 
of stellar mass and the duration of the SMG phase. 

The median value of the SFR per unit stellar mass 
(SSFR = SFRir/A4, Column 7 of Table ESj of ~ 
1.8Gyr~^ is within the range for other high-z star- 
for ming samples (compare wi th Figures 2 and 4 of 
iCastro Ceron et all . l2006l l2009l respectively). This indi- 
cates that SMGs are forming stars intensely. 

SSFRs are compared with (the inverse of) the age of 
the Universe in Figure H] The SMGs close to the solid line 
could have formed their stellar populations at the present 
rate within the age of the Universe. However, the SMGs 
close to, or above the dashed line could have formed their 
stars at the present rate within less than 10% of the age 
of the Universe, i.e., within < 300 Myr at z = 2. These 
galaxies are experiencing a powerful starburst episode. 

At the extreme there are three high-z SMCfo with very 
high SSFRs > 10 Gyr-i (Column 7 of Table 1X31). They are 
all hot {Td > 60 K, Column 13) and formed the majority 
of their stellar populations during the ongoing starburst 



SMMJ131201.17±424208.1, SMMJ141802.87±523011.1, 
SMMJ221806.77±001245.7 plus a low-mass, low-2 case, 
SMMJ030238.62±001106.3 



8 



Michalowski et al.: Evolution of submillimeter galaxies 



(Mburst/M* > 60%, Column 9). Therefore they are likely 
the most powerful cases of SMGs formed in major mergers 
of galaxies with huge gas reservoirs that were subsequently 
converted into stars. 

Our median SSFR at z > 1.7 (1.8 3 Gyr~ ^) is a factor 
of ~ 2 lower than that of iDunne et"all (|2009l . 3-4.5 Gyr"^; 
see their Figure 12b) for 10" < < 10^^ Mq galaxies at 
these redshifts. This difference can be ex plained if the radio 
luminosities (used bv lDunne et al.l . [20091 to estimate SFRs) 
are boosted by AGN activity more than the IR luminosi- 
ties used here. Indeed, if we use SFRradio instead of SFRir 
to calculate SSFRs the median for the SMGs at z > 1.7 
increases to 3.20 Gyr~^ (see Section [5.4.21 for discussion of 
AGN contamination in our sample). 

In order to assess the accuracy of SFR estimates based 
on radio emission (independent of SED modeling) we com- 
pared the ratio of SFRj-adio/SFRiR. Its median value is 
equal to 1.3. Hence, assuming that IR emission is a good 
proxy for SFR, then radio estimates suffer from a ~ 30% 
systematic error. This is illustrated on Figure [5] where the 
dashed line denotes the relation between IR and radio lu- 
minosities required to make SFRir = SFRj-adio- Indeed the 
radio luminosity gives systematically higher SFRs for SMGs 
(most of the points are above the line). This can be caused 
by a significant AGN contamination boosting radio flux (see 
Section [5.4.2p . or a strong bias favouring radio-bright galax- 
ies, because those non-detected at radio do not enter our 
sample (Section[2]). Alternatively, it coul d be that fo r lumi- 
nous galaxies either the IR conversion ofE enmcu t3 (IT998h 
should be scaled u p by a factor of 1.3, or the radio conver- 
sion of mi ()2003[ ) scaled down. 

5.2.2. Stellar masses 

SMGs h aving stellar masses of lO^-lO^^ Mq (Column 8 
of Table IA.3I and Figure [5]) are among the most massive 
galaxies in the Universe, regardless of redsh i ft (compare 
with Figures 2 and 4 of ICastro Ceron et all . l2006l l2009l 
respectively) . This property makes them natural candidates 
for the progenitors of the present-day ellipticals. 

The relatively tight range of stellar masses is likely not 
a result of sensitivity limits at optical and near-IR. This is 
because i) galaxies with stellar mass as low as ^ 10^ Mq 
would have been d etected in de e p Spi tzer imaging at red- 



shifts z ~ 2 (e.g. iReddv et all I2006D ii) our samp le ac 



sam- 



counts for 50% of the parent IChapman et all (pOa 
pie (and only 30% of the parent sample may have different 
properties than our sample, see Section [5]) , so it is unlikely 
that we miss only the low-mass objects. Therefore, high M* 
seems to be an intrinsic property of submillimeter-selected 
galaxies. Mergers of less massive galaxies could not result in 
a powerful starburst giving rise to detectable submillimeter 
emission (see also Dave et al., 2009). 

Only a minor part (median ^ 8%, Column 9 of 
Table IA.3I and Figure [2]) of the stellar populations present 
in SMGs has been formed during the ongoing starburst 
episodes. Hence, even though SMGs probably evolve into 
ellipticals, the majority of the stellar mass in such ellipticals 
had been created before the submillimeter-bright phase. 

This could mean that the current SFRs and stellar 
masses of SMGs are only loosely connected and indeed this 
manifests itself in a very high spread (around two orders of 
magnitude) in SSFRs in our sample even though the stellar 



mass range is relatively tight: ^ lO-'^^-lO^^ Mq (Figure [2). 
Th is behaviour is unusual compared to other galaxies (see 
Ca stro Ceron et al .. 2006, 2009). 

However we note that the low stellar masses created 
in the ongoing starburst may partially be an effect of the 
assumed starburst ages of 50 Myr. If a starburst dura- 
tion of 100- 200 Myr were adopted ( ISmail et al.', '200^ 
iBorvs et all . |2005; Hainhnd . l2008t iTacconi et al. . ,2008) the 
resulting Mburst could be higher by a factor of ^ 2-4. 

The mass- to-light ratios, M^/Lk, of SMGs (Column 10 
of Table IA.3I and Figure ^ are typical for massive galax- 
ies. Specifically, the median (0.68 Mq L|^^) is sim ilar to 
the values for M* > 10" Mq galaxies ()Drorv et all . |2004 
their Table 1) and t o simulated massive galaxies at ~ 1 
(jCourtv et al.l l2007l their Figure 4). 



5.2.3. Dust properties 

Our fits suggest that SMGs are moderately dust-obscured 
with a median Ay ^ 2 mag (Column 14 of Table IA.3|) . 
Our estimates are consistent withi n l-2 i7 with the 
mean/median values obtained by Small et al. ([2004 , 1.70- 
2.44). ISwinbank et al (2004. 3.0 ± l.O'lT lBorvs et al.l (llOOl 
1.7 ± 0.21 and iHainline (.20081 1.7 ± 0.1) based on near- 
IR data. For individual SMGs we obtained systematically 
larger ext inctio n (median difference of ~ 0.3 mag) than 
Hainling (I2008D . The difference may be accounted for if 
there is significant extinction even in Spitzer IRAC data. 

The dust density of SMGs at z < 0.5 (Column 9 of 
Table [IJ is approximately 3% of the total local (0.013 < 
z < 0.18 ) dust budget of log pdust — 5.57t°:l?MQMpc " 



given by iDriver et all (|2007l) based on an assumed dust- 
to-light ratio. Therefore SMGs contribute very little to the 
dust budget at low redshifts. 

In our sample of SMGs pdust does not change signif- 
icantly from z ~ 3.6 to z ~ 0.5. We do not detect any 
evolution of dust mass in SMGs across the entire redshift 
range (Figure [2]). A constant dust mass density a cross red- 
shifts 0-3.5 was also found bv lPascale et al.l (j2009) based on 
a stacking analysis at submillimeter wavelengths of galaxies 
selected at 24 /xm. 

The question is what happened to the dust produced 
in SMGs. If they evolve into dust-poor ellipticals, then the 
dust is not simply stored in their end-products (as is proba- 
bly the case for stellar masses). It is therefore plausible that 
dust is either blown away (by stellar and/or AGN winds) or 
absorbed in star formation, or destroyed during subsequent 
evolution after the SMG event. 



5.2.4. Comparison with GRB hosts 
In I Michalowski et al.l 



we presented a hypothesis that 
gamma-ray burst (GRB) host galaxies may constitute a 
subsample of hotter/less luminous counterparts of SMGs. 
Indeed, the UV-to-IR SEDs of three z ~ 2-3 SMGE are 
consistent with z I submillimeter /radio bright GRB 
hosts (dashed lines on Figure EHI from [Michalowski et all 
I2008D . but 1.2-3.9 times more luminous. These three SMGs 
are similar to GRB hosts with respect to their hot dust 
temperatures (~ 40-60 K), high SSFRs (> 2 Gyr"\ high 



* SMMJ141750.50-H523101.0, 
SMMJ163627.94+405811.2 



SMMJ141802. 87+523011.1, 
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Fig. 4 Top: Cosmic star formation density. The SMGs' contribution rises with redshift from ^ 9% to 20% fSection l5.3.1l 
and Table [T]). Filled Squares: data for SMGs at z > 0.5 in four bins (Table [T] and Sectioning)- Small Squares: data 
for SMGs at z < 0.5. Thick black arrows: the SMG data without incomp leteness correction (fac tor o f 3.5, Sectio n l3T2l) . 
Black crosses and diamonds: star formation density of SMGs determined bv lChapman "eTall (pOOSi) and lWall et all (|2008D . 
respectively. Colored points with error bars: determination of the cosmic value based on different estimates - ultraviolet 
{violet), emission lines: [O 2], [O 3], Ha, H/3 (green), mid-IR {light blue), submillimeter {black), radio {red). X-ray {yellow). 
Extinction correction and, in many cases, incompleteness correction have been applied by the authors. Arrows: lower 
limits. 

Bottom: Cosmic stellar mass density. The SMGs' contribution rises with redshift from ^ 5% to ^ 50% f Section 15.3.21 
and Table [1]). Red points with error bars: determination of the cosmic value from literature. The data and the references 
are listed in Tables IA.4I and IA.5I in appendix. 
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fraction of stellar mass formed in the ongoing starburst 
(> 10%) and blue optical colors. 

If larger samples of GRB hosts shows a similar tendency 
that their brightest members overlap with the hotter sub- 
sample of SMGs, then GRB events will provide an effective 
way of selecting hot SMGs, otherwise difficult to localize. 

5.3. Contribution to stellar mass assembly 

5.3.1. Star formation rate volume density 

SFR densities of SMGs were calculated as described in 
Section 13.21 In order to assess the accuracy of our sim- 
plified method of dividing the sum of the SFRs of the 
detected SMGs by the total survey volume, we compare 
our estimates with those resulting from detailed calcula- 
tio n of the volume contribution of individual SMGs done 
by iChapman et aL (f2005„ based o n the same sample as 



we analvse) and IWall et al] (|2008l based on 35 SMGs in 
GOODS-N field of which 17 have spectroscopic redshifts). 
The comparison is shown in Figure 01 Our results in two 
high-redshift bins (z > 2) corrected for incompleteness 
(Section 13. 2 [) are cons istent with that of IChapman et al.l 
(|2005D and IWall et all (I2008D. At lower redshifts we find 
values similar to 'Chapman et al.' (2005), but an order of 
magnitude lower than Wall et al. (2008). Therefore we con- 
clude that i) our method to calculate volumes is accurate, 
since it gives consistent results with other estimates; and 
ii) our sample is incomplete in the three low-redshift bins 
as anticipated in Section |21 

From Figure HI (and Columns 4 and 5 of Table |T]) it 
is apparent that a psFR of SMGs starts to decline (with 
cosmic time) earlier ( about z ~ 2) than that of other 
galaxies (z ^ 1). More quantitatively, SMGs harbour ^ 
20% of the cosmic psFR at z 2-3.6 (Column 5), but their 
contribution drops to ~ 9% at 0.5 < z < 1.4. It is likely that 
at lower redshifts, due to the decreased ra. te of mergers (e.g. 
iRawat et all . 120081: Ide Ravel etaP . I2009D . there are fewer 
galaxies left that can still sustain high SFRs to be detected 
at submillimeter wavelengths. However, part of the decrease 
of SMG psFR can be explained by the "redshift desert", 
which makes it difficult to measure redshifts of z ~ 1.2-1.8 
SMGs (see Section 

A high value of psFR of SMGs at z ^ 2-3 and 
the subsequent decline are consistent with the hypothesis 
that the SMG population is a manifestation of powerful 
starburst episodes evolving into the present-day ellipticals 
(as discussed in Section I5.2.2p . In this scenario galax- 
ies detected in the submillimeter at high-z do not en- 
ter the sample of SMGs at low-z because they have al- 
ready evolved into passive galaxies. It has indeed been 
found that elliptica l s cont a in old stars formed a t z ^ 
1.5-4 (iDaddi et al.'. 2000; van Dokkum fc Franxl . l200lt 
Ivan de Ven et all . ,2003 ). The evolution of SMGs into el- 
lipticals has also been clai med by sever a l auth ors based on 
their luminosity function (jSmail et al.l . l2004[ ). huge lumi- 
nositics (Ealcs ct al., 199^ and gas reservoirs (jSmail ct aL, 
2002; G reve et al.l. I2005D. strong clustering (|lvison et al. . 
2000; .Almaini et al.L 120031) . space density and morphol- 
Qgv |Barger et all. 119991: iLillv et al. l.ll999t(Tre ntham et all . 
11999*: 'Swinbank et al.l. l2006D and evolutionarv SED models 
(T akagi ct al., 200 4). 

iKnudsen et al.l (|2008b( ) analysed number counts of 
SMGs fainter than the SCUBA confusion limit, using those 



behind clusters of galaxies magnified by lensing. They con- 
cluded that the integrated light produced by the SMGs 
brighter than 0.1 mjy (i.e. LIRGs and ULIRGs with 
roughly Ljr > 8 x 10^° Lq and SFR > 15 Mq yr^^) is 
comparable to t he extragala ctic b ackground light (EB L) at 
850 Atm (see also lBlain et aLl . ll999t[Cowie et al.l . l2002D . This 
means that these galaxi es host the majority of the cosmic 
obscured star formation. IKnudsen et al.l (|2008b| ) also found 
that sources brighter than 2.5 mJy (roughly the limit of the 
survey considered here) contribute ^ 25% to the to EBL 
at 850 nm. (see also Hughes et al.l ll998l : iBarger et al.l . 119991 : 

Together with our 
80% (4x20%) of the 



IWang et all |2004 ICoppin et al 
results this implies that as much as 
cosmic star formation at z ~ 2-3.6 reside in SMGs brighter 
than 0.1 mJy. This is only true if the faint (< 2 mJy) SMGs 
have similar dust temperatures to the brighter ones. If they 
are colder (hotter) their submillimeter fluxes corresponds 
to lower (higher) SFRs (because it is calibrated to total 
IR emission) and therefore the total SMG population con- 
tribute less (more) than 80% to the cosmic psFR- This 
picture is however complicated, because based on stacking 
analysis it has been claimed that the distri bution of the 
faint SMGs peaks at lower redshifts (z < 1.5: IWang et al.l 
[20061: ISCTieant et al.L 120081) . 

Our overall conclusion is that the SMG population plays 
a significant role at redshifts z ~ 2-4, namely sources 
brighter than - 3 (0.1) mJy at 850 /im host 20% (80%) 
of cosmic star formation. Their contribution can however 
be lower in reality if very small (but numerous) galaxies 
are missed in all high-z fiux-limited galaxy surveys. In such 
a case the total SFR density (color points on Figure SI) 
would be underestimated. To solve this issue much deeper 
surveys at high-z are necessary, either blank-field or for 
well-selected dwarf galaxy samples (e.g., GRB hosts or Lya 
emitters). 

IZheng et al.l (|2007l) estimated psFR at z 0.9 for mas- 
sive galaxies (M* > 10^^ Mq) down to i? < 24 mag 
40% of SMGs satisfy the latter criterion) equal to 



(only 

0.0052j:°;|^°20 yj.-i Mpc"^. This value is only a factor 
of 2 lower than our estimate for the SMGs at 0.5 < z < 1.4 
(Table [H). Therefore, although SMGs do not host a ma- 
jor fraction of the cosmic SFR at these redshifts, they con- 
tribute significantly (0.0102/(0.0052 x 0.6 + 0.0102) - 66%) 
to the SFR budget of massive galaxies. 



5.3.2. Stellar mass volume density 

Stellar mass densities of SMGs were calculated as described 
in Section [XH Figure H and Table [H (Columns 6 and 7) 
show that at z ^ 2-3.6 a significant part (~ 30-50%^ of 
the cosmic stellar mass had been formed in the progenitors 
of SMGs. At lower redshifts p* of SMGs (and hence their 
contribution to the cosmic /?*) drops, likely because the 
majority of SMGs at higher redshifts had already evolved 
into passive galaxies at z ~ 1.5, and so dropped out of 
our submillimeter-selected sample. Moreover the sample is 
incomplete at z ^ 1.2-1.8 due to the "redshift desert" (see 
Section 121). This brings down the densities of SMGs in the 
low-z bins. 

Since most of the stellar mass of SMGs has not been 
formed in the ongoing starburst fSection I5.2.2p . their 
reflects the integrated contribution of SMGs to the cosmic 
PSFR- Therefore the relatively high contribution of SMGs 
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to the cosmic in the last redshift bin 31%, Column 7 
of Table [T]) means that SMGs play a non-negligible role 
in the cosmic stellar assembly even at z > 3.6. This can 
be tested by analysis of a sa mple of z > 4 SMGs in 
a defined survey sky area (e.g. I Michalowski et all l2010t 
lYounger et al.l . l2009al , note that these results are likely af- 
fected by cosmic variance). It has b een confirmed th at such 
distan t SMG s exist (Capaket al. . [ 2008: Knudsen et al 
2008al \20ldi . Schinn erer et al.L l2008n ( 
Daddi et al.f . l2009bi.al) . 



Coppin et all . 12009 : 



5.4. Source of emission 
5.4.1. IR-radio correlation 

With our full SEP modelling of 76 SMGs we confirm the 
results of iHainlind (|2008[ ) on the correlation between IR 
and radio luminosities. Figure [S] shows that SMGs follow a 
linear IR-radio correlation. The two outliers (with q ^ 1.3, 
see Section 15.4.21) are probably caused by AGN activity 
contributing significantly to radio luminosities. A linear fit 
gives: 

log(L,i.4GHz/ioHz-i) = (0.95 ± 0.07) log(iiR/L0) 

-(14.3 ±0.8). ^ ' 

The slope is consistent (within errors) with unity, suggestive 
of the linear relation between L^, 1.4 ghz and Lir at the high- 
end (LiR > 10^^ Lq) of the galaxy lumin osity fu nction (a 
similar value of 1.064±0.025 was found bv lHain hnc. 2008). 

The IR-radio correlation is usually quantified by the 
ratio of IR and radio luminosities, q (see Section [5TT|) . The 
mean q for SMGs (2.32 ±0.04, scatter: 0.34) is significantly 
lower than that of local star-forming galaxies (2.64 with a 
scatter of 0.26; feell, 20 0j). Similar o ffsets were reported 



bvlKovacs et al 



. . l2006D . lMurphv et al.l ( 20091 and Murpl^ 

()2009() based on smaller samples of SMGs. We conclude 
that at z > 1.4 SMGs have radio luminosities on average a 
factor of ^ 2.1 larger (Aq ^ —0.32) than what would result 
from the local relation. The difference is significant at the 
level of 4-5(7 and can be explained in three ways. 

Radio-loud AGNs have on ave rage low q values (see e.g . 
iMiUer fc OwctI . [200TI: lYun et all . [200lt lYang et al.L [200I . 
If > 50% of the radio emission of SMGs is powered by 
AGNs, then the radio luminosities of SMGs higher by 
a factor of ^ 2.1 can be accounted for. However, there 
are indications that SMGs are starburst-dominated (see 
Section [5. 4. 2p . so we deem this explanation less likely. 

Another explanation is that the radio excess is a re- 
sult of the bias against radio-faint sources in our sam- 
ple (see Section [2]) . This can be tested when a sample of 
SMGs with localizations (and hence redshifts) in dependent 
of radio detect i ons is available (e.g. Daddi et a l.. 2009b3; 
iKnudsen et all . 120101: IWeifi et al.l . l2009all . 

The third possibility is that some properties influ- 
encing the IR or radio emission are intrinsically d iffer- 
ent fo r SMGs and local galaxies. The sample of I Belli 
(|2003f) includes local normal, star-forming spiral and ir- 
regular galaxies, blue compact dwarfs, starburst galax- 
ies and ULIRGs. Therefore the difference in the prop- 
erties between this sample and such extreme galaxies 
as SMGs i s expe cted. Such explanation w as offered by 
iLacki et al.l (|2009f ) and iLacki fc ThompsonI (p009 ). Their 
numerical modelling showed that cosmic-ray electrons 
in "puffy starbursts" (vertically and radially extended 



galaxies with vertical scale heights ~ 1 kpc) experi- 
ence weaker bremsstrahlung and ionization losses result- 
ing in stronger radio emission. Indeed, there are indica- 
tions that SMGs are ex t ended on vertical scales of ~ 1 
kpc (ILacki fc ThompsonI. 120091: iTacconi et al.l [20061 120081: 



iGenzel et all . 120081: lYounger et al.l 120081: iLaw et "311120091 ) . 

so we find this explanation probable. 

The systematic uncertainties in the determination of 
LiR (factor of < 2, Section Uj) may in principle also ex- 
plain the offset. However, we find this unlikely because 
similar offsets w ere found by other authors using different 
fitting inethod s (|Kovacs et al.l 120061: iMurphv et all . 120091: 
lMurphvll2009l) . 

The q values for SMGs are shown in Figure |S] as a func- 
tion of redshift. We do not detect any significant evolution 
across the redshift range 1.4-3.6. The only sign of evolution 
is that the mean q in the low-redshift bin (0.5 < z < 1.4) is 
above the value found at higher redshifts (~ 4tT). This can 
be explained either by the contribution of reproc essed emis- 
sion from low-mass stars (cirrus emission, e.g. lYun et all 
,2001, and references therein) to the IR, or by the fact that 
at low redshifts SMGs are more similar to other local galax- 
ies and do not exhibit large vertical scale heights charac- 
teristic for "puffy starbursts" (see above). 

It is important to note that the derived linear IR-radio 
correlation for SMGs is not a consequence of the use of the 
SED templates (which were tuned to fulfill this correlation 
locally), because the radio luminosities used here were de- 
rived based on the observational data only, independent of 
the SED modeling. 

5.4.2. AGN activity 



As discussed in Section r5. 4. 11 AGN activity could explain 
low q values of SMGs. This is at least true for the two 
SMGs with lowe s t jj, spectroscopically classified as AGN 
(|Chapman et al.l . l2005D . 

In the SEDs of SMGs there are clear signs that some of 
them host AGNs (though, not necessarily a bolometrically 
dominant ones). Radio datapoints are higher than model 
predictions by more than 3cr in 36% (27/76) of SMGs, 
whereas they are lower than models only for 8% (6/76). 
This may hint at an AGN contribution in these galaxies. 
However, 4 out of 5 X-ray identified starbursts (Column 16 
of Table ETSl) also exhibit radio excess, so we find other ex- 
planations of radio excess presented in Section 15.4.11 more 
reliable. 

Another indication of an AGN contribution is that 18% 
(14/76) of SMGs show a mid-IR power-law AGN feature in- 
compatible with our starburst models (see Figure lA.ll and 
Column 16 of Table IX3|) . However, rest-frame 2-5 /im ex- 
cess was also interpreted as a tracer of recent star formation 
(Mcntuch ct al., 2009). 

Finally, three SMG£EI have exceptionally high SFRuv 
(> 500 M0 yr-\ Column 4 of Table |X3]). Strikingly, all 
of them were fitted with non-starburst models (Mburst = 
0, Column 9), so modeling is consistent with these high 
SFRs being continuous (the same is true for three other 
non-starburst SMGs with high SFRir). Such a scenario is 



^ SMMJ131215. 27+423900.9, SMMJ141813. 54-1-522923.4 
^° SMMJ123716.01+620323.3, SMMJ131215. 27+423900.9, 



SMMJ131222. 35+423814.1 
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Fig. 5 Radio luminosity density as a function of infrared (8-1000 /im) luminosity of SMGs showing a linear relation, though 
with a normalization offset from the local relation by a factor of ~ 2.1 towards higher radio luminosities fSection l5.4.1[) . 
Circles: values for individu al SMGs color-coded by redshift. Solid line: linear fit to the data (eq.[Il). Do t ted lin e: the mean 
local rela tion (iBelll. 20031) . Shaded area: its scatter. Dashed line: the track where SFRir (lKennictit3 . 1 19981 ) is equal to 
SFRradio (|Beli 12003 ). The strong outliers (above the line) at high- luminosity end are probably caused by AGN activity 
increasing radio luminosities. 



unlikely, so this hints at an AGN contribution to the UV/IR 
emission. 

However, the fact that we obtained reasonable SED fits 
for most of the SMGs using purely star-forming models 
(Figure ET} hints at the conclusion that AGN activity is 
not dominant in our sample. 

We investigated the issue of AGN activity further by 
analysing the average q values of the follo wing subsamples 
(see also Figure[6]): X-ray identified (Alexander et al.l . 120051 ) 
AGNs: 2.32± 0.06 and starbursts: 2. 12±0.18; optically iden- 
tified AGNs (|Chapman et al.l . l2005D : 2.27 ±0.09; and mid- 
IR identified AGNs (see above): 2.36±0.12. All subsamples 
are consistent with the value derived for the entire sam - 
ple (2.32) Hence, we confirm the finding of lHainlin^ (|2008[ ) 
that even the AGN-classified SMGs follow a linear IR-radio 
correlation. This means that even if an AGN is present it 
does not contribute to the emission of an SMG significantly 
(with the exception of the two g ~ 1.3 sources). 

This is in line with i) the X-ray studies of SMGs indi- 
cating that the contribution of A GN activity to their IR 
emission is only ^ 8% on average ([Alexander et al.L l2005[) ; 
m) mid-IR colors of SMGs indicating that AGNs dominate 
the emission at these wavelengths only in 13-19% cases 
(jHainline et al.L l2009l) : Hi) mid-IR spectroscopy of SMGs 



reveal ing only w eak AGN-like contin ua (iValiante et al 



Pope et all. 2008; Mcncndez-Dcl mestre et al.l . 12007 



Murphv et all . 120091: iWatabe et all . l2009t ): iv) 



near- 

IR spectroscopy rev ealing that starbursts dominate the 
emission of SMGs (ISwinbank etaD . I2004D . Moreover, 
Ide Vries et al.l (|2007l ) found that star formation processes 
(if present) account for at least 75% of the radio luminosi- 
ties of optically-selected AGNs. 

Therefore we conclude that AGNs are present in a sig- 
nificant fraction of SMGs, hut their contribution to the IR 
emission is at most minor. 



5.5. Comparison of our results with the literature 

For the sample of SMGs discussed in this paper there are 
previous estimates of some of their properties. In this sec- 
tion we compare them with our results. 

IChapman et al.l (|2005l ) derived Ljr and Td based only 
on the 850 /Ltm and 1.4 GHz data. There is no systematic 
difference between the determinations of Td (our median 
of 38.7 K, theirs: 38.3 K). The mean difference between 
individual datapoints is 4 K (■^ 10%). However, our val- 
ues for LiR are systematically lower than theirs (the me- 
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offset from the local relation by a factor of ^ 2.1 {Aq ~ —0.32) towards higher radio luminosities (Section l5.4.1l) . Circles: 
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dian ratio of individual datapoints is 1.7). We find our val- 
ues more reliable since they are based on data spanning a 
wider wavelength range. Overestimation of L ir when using 
only 850 /im and 1.4 GHz was also noticed bv lKovacs et al.l 
12006.) and Pope et all (|2006D . 

Kovacs et al .1 (|2006D investigated a subsample observed 
at 350 /im. Their median dust mass (9.04 logM©) and q 
value (2.20) are consistent with our estimates (9.01 and 



The median difference between individ- 
30% for dust masses and ^ 13% for 



2.35, respectively) 
ual datapoints is 

q- 

The median stellar mass for a subsample of 13 SMGs 
investigated by .Borys et al. (2005, 11.51 log A/q) is close 
to our value (11.70). However, estimates of iHainlind (j2008l 
median 10.82 logM©) for 64 SM Gs are a f a ctor o f ~ 5.6 
smaller than our values (11.57). iHainUnd EoM) postu- 
lated that the disc repancy between her results and those of 
Hor^sliaD (|2Q0l) arose from a combination of systematic 
differences between the applied SED models and a higher 
AGN contribution in the K-hand (used by iBorvs et al.l . 



l2005f l with respect to the iJ-band. Our estimates are based 
on all the available photometric data, and so we find the 
former explanation more likely. In particular, the differ- 
ences in the applied stellar population models and their 
ages may explain the discrepancy. 

6. Conclusions 

We have investigated the UV-to-radio SEDs of 76 SMGs 
(£"850 ^ 3 mJy) with spectroscopic redshifts (0.080-3.623). 
For the first time the properties of such a significant sample 
has been derived consistently using all available data. The 
resulting SFRs (median 713 M© yr~^) and stellar masses 
(11.57 logAf©) are among the highest in the Universe. 

Such high stellar masses, already present at redshifts 
^ 2-3, require that SMGs experienced either at least two 
star burst episodes, or a merger of several smaller galax- 
ies. Our modeling suggests that only a minor fraction (8%) 
of their stellar populations was formed during the ongoing 
starburst episodes. This is supported by the fact that the 
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SFRs and of SMGs are basically disconnected, i.e. we 
observe two orders of magnitude spread in SSFRs whereas 
the range of M* is relatively narrow: 10^^-10^^ Mq. We 
concluded that dust is blown away or destroyed during the 
evolution of SMGs, since it is not stored in the likely end- 
products of SMGs, elliptical galaxies. 

Indeed, the high stellar masses and the evolution of the 
SFR and stellar mass densities of SMGs are consistent with 
a scenario in which SMGs are progenitors of present-day 
ellipticals. 

We found that SMGs contribute significantly to the cos- 
mic SFR, psFR (~ 20%) and stellar mass, (30-50%) 
densities at z ~ 2-4. If we consider submillimeter sources 
down to 0.1 mJy the contribution to psFR rises to ^ 80%. 

Our analysis suggests that a linear IR-radio correlation 
holds for SMGs at least up to a redshift of 3.6, but they 
are '^2.1 times brighter at radio wavelengths than what 
would result from the local correlation. 
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Table A.l Photometry detections of SMGs 









flux 


error 




SMG 


z 


(/xm) 




(^Jy) 


Reference 


SMMJ030226. 17+000624. 5 


0.080 


0.365 


33.419 


0.306 


Clements et al. (2004) 


SMMJ030226. 17+000624. 5 


0.080 


0.428 


83.176 


7.319 


Chapman et al. (2005) 


SMMJ030226. 17+000624. 5 


0.080 


0.440 


83.946 


0.232 


Clements et al. (2004) 


SMMJ030226. 17+000624. 5 


0.080 


0.550 


181.970 


0.335 


Clements et ah (2004) 


SMMJ030226. 17+000624. 5 


0.080 


0.656 


275.423 


24.234 


Chaoman et al. (2005) 



Note. — The second reference in the first entry for each SMG indicates that the re exists a mid-IR Spitzer/IRS spectrum of 
this object (jValiante et all . l2007l : iPooe et"all . l200^ : iMenendez-Dehnestre et al.l . 12007. 2009). This table is available in its entirety 
in a machine-readable form in the electronic edition of the Journal. A portion is shown here for guidance regarding its form and 
content. 



Refe rences. — llvison et all (I20021. 12 005): Chapma n et all ll20 03b. 2005): C apak et all ll200l ): IClements et al. ' (2004): Egami et al 
2003l : |Greve et al.l ll2004l ') : ISmail et al.. (,2004.' ): | Fomalorit et all (|200& '): .Kovacs et al.l (f2006l '): iLaurent et al.l (j2006l : .Tacconi et al 
200^ : lPope et all ((200^ : iHuvnh et aLl (pOOTl ) : iHaiiilind (|2008l ) 



Table A. 2 Photometry upper limits of SMGs 







Aobs 


flux 


error 




SMG 


z 




(MJy) 


(MJy) 


Reference 


SMMJ030227.73+000653.5 


1.408 


70.000 


13600.000 


0.000 


Mainline (2008) 


SMMJ030231. 81+001031. 3 


1.316 


0.365 


0.059 


0.000 


Clements et al. (2004) 


SMMJ030231. 81+001031. 3 


1.316 


0.428 


0.100 


0.000 


Chapman et al. (2005) 


SMMJ030231. 81+001031. 3 


1.316 


0.440 


0.102 


0.000 


Clements ct al. (2004) 


SMMJ030231. 81+001031. 3 


1.316 


0.656 


0.331 


0.000 


Chapman et al. (2005) 



Note. — When the error is equal to zero, the flux column denotes 3cr upper limit. Otherwise — formal flux at the position of an 
SMG. This table is available in its entirety in a machine-readable form in the electronic edition of the Journal. A portion is shown 
here for guidance regarding its form and content. 



Refe rences, —llvison et al . ( 2002, 2005); Chapma n eraP (l2 003b'.'2005'):'C apak et al ' (2004"): 'Cle ments et al .' (2004"): 'Egam i et al 



2004 ) 



2006D 



Greve et al.l (120041): iSmail et al (.2004) : Fomalont et al (200& ): .Kovacs et al.1 (2006 ): .Laurent et al.l (|2006. ): Tacconi et al 
Pope et al.1 (|2006D : iHuvnh et all (120071 ): iHainlind (|2008l ) 



Table A. 3. Properties of SMGs derived from the SED modeling 



SMG 


z 

(2) 


SED 


Sll! [M 
UV 
(4) 


vr M 
ii\ 

(0) 


— radio 

f(\\ 
(o) 


SSFR 
(Gyr-i) 
(7) 


logM, Mwst/M. 

(Mo) (%) 
(8) (9) 


MJLk 

(Mq/Lq) 

(10) 




log LiR 

(Le) 

(12) 


Td 

(K) 
(13) 


Av 

(mag) 
(14) 


9 
(15) 


AGN*^ 

(W) 


Diviivi.J UoUzz / . / v> + UUUOOO.O 


1 


408 


825 


7 


91 


713 


1 190 


4 


81 


11 


17 


25 


5 





42 


8.54 


12 


62 


52 


6 


1 


81 


2 


29 


rad 


DiVliVi.JUoUzoi.oi+UUiUoi.o 


1 


316 


433 





68 


224 


211 


7 


81 


10 


46 


71 


1 


1 


36 


9.10 


12 


11 


25 


6 


4 


38 


2 


54 




oiVliVi.J UoUzoD. iO+UUUoi / . i 


2 


435 


91 


4 


18 


1152 


810 


3 


33 


11 


54 


1 


1 





55 


8.62 


12 


83 


41 


3 





27 


2 


67 


mllx,rad 


OiViiVi.J iUOiOi.Dy+O / ZDOD.U 


1 


147 


134 





30 


248 


457 





92 


11 


43 


1 


5 





65 


8.97 


12 


16 


34 


4 


4 


00 


2 


25 


rad 


oMMJ iUoioo.4 f +0 r zoiz. r 


2 


686 


837 


12 


95 


589 


1110 


1 


91 


11 


49 


12 


5 





66 


8.90 


12 


54 


38 


7 


2 


29 


2 


24 


mIR. 


oMiVlJ iUOloo.Uz+O* ioUU.z 


2 


239 


314 


16 


77 


373 


1561 





44 


11 


93 





6 





60 


9.15 


12 


34 


33 


2 


1 


12 


1 


89 


rad 


QAyTAyr Ti nc;onn oo_i_c;?oyion o 

oiViiViJ lUOZUU.ZZ-|-0 1 





689 


1 HQ 

loo 





26 


1 1 ft 
lies 


oy 





74 


11 


20 


4 


5 


1 


56 


ft TO 


11 


84 


33 


2 


2 


58 


2 


82 




oiVilVi J iUozUi .Z0-|-O / ZHl^O. 1 


2 


148 


yoo 


2 


58 


798 


1044 


3 


40 


11 


37 


18 


1 





95 


9.01 


12 


67 


45 


1 


2 


91 


2 


40 


inlK. 


QlV'f AT T1 fT^OriT /I 0_l_ 1^71 Qfl/I n 


2 


689 


2 195 


3 


81 


1559 


6679 


1 


93 


11 


91 


12 


5 


1 


18 


8.80 


12 


96 


45 


1 


2 


85 


1 


88 


rad 


Q A 'f A T T 1 n ^OO 70_i_f^71 onK /I 
OiViiVi.J iUOZZO. i\)-\-o I ±yUD.4: 


2 


372 


217 


18 


89 


329 


2310 





75 


11 


64 











49 


8.96 


12 


28 


33 


2 


1 


16 


1 


67 


rad 


OiViiVi.J iUozz i .oo-|-o t ZOiZ.4: 


2 


142 


410 


4 


73 


376 


564 





57 


11 


82 


1 


9 





78 


8. 79 


12 


34 


38 


7 


1 


71 


2 


34 




DiViiVi J iUOZZ ( . ( (-|-0(ZZ±O.Z 


1 


956 


515 


1 


55 


A AT 


A T'i 

4 r o 





89 


11 


70 


4 


5 


1 


16 


9.01 


12 


42 


33 


2 


2 


17 


2 


49 




oiViiVij iuozou. ( o-|-o ( zzuy .O 


2 


611 


1255 


45 


92 


975 


iy4z 





72 


12 


13 


3 


8 





68 


9.18 


12 


75 


38 


7 


1 


51 


2 


22 




OiViiViJ iUOZoo. ly + O / iDOi. i 


1 


852 


919 


21 


87 


659 


734 


4 


54 


11 


16 


28 


8 





37 


8.85 


12 


58 


38 


7 


2 


13 


2 


47 




oMiVlJ iUozoo.oU+O / z4oo.o 


3 


036 


1558 


7 


40 


1169 


1921 





69 


12 


23 


3 


9 





85 


9.04 


12 


83 


45 


1 


1 


86 


2 


30 


spec 


OiViiViJ iZo04y.44-|-DZiOOO.O 


2 


203 


237 


29 


23 


335 


1144 





70 


11 


68 











51 


9.18 


12 


29 


33 


2 


1 


11 


1 


98 


X,rad 


oiViiVi J izoooo.zD+Dzioo / . / 


2 


098 


1460 


6 


05 


1223 


802 


5 


73 


11 


33 


31 


8 





42 


8.85 


12 


85 


52 


6 


2 


99 


2 


70 




OiViiViJ izoOOO. i4-|-DZUyUi . / 


1 


875 


332 


11 


94 


358 


2252 





61 


11 


76 


1 








61 


8.98 


12 


32 


33 


2 


1 


41 


1 


72 


A.,mlll,raa 


OiViiViJ izoDUU. iU+ozUzoo.o 


2 


710 


3733 


6 


49 


3066 


6407 


8 


92 


11 


54 


51 


6 





29 


8.52 


13 


25 


71 


5 


3 


64 


2 


20 


rad 


OiViiViJ izoDUU. io-|-OziU4 / .z 


1 


994 


1449 


3 


04 


1102 


1602 


2 


53 


11 


64 


15 


2 





89 


9.03 


12 


81 


38 


7 


2 


80 


2 


35 


SB, rad 


OiViiViJ iZOuUD. ( Z-|-DZiOOU. / 


2 


416 




30 


23 


354 


Ar^A 
404 


1 


27 


11 


45 


6 


9 





65 


8.92 


12 


32 


33 


2 


1 


15 


2 


41 


X,spec,mIR 


OiViiViJ lZoDUD.oO-|-OZlUZ1.4 


2 


509 


1 OQ1 


12 


59 


1 1 K'i 

1 loo 


iooi 


1 


17 


11 


99 


5 


6 





79 


y.o4 


12 


83 


33 


2 


1 


58 


2 


39 


A.,raci 


OiViiViJ IZoDlD. iO-|-DZ10iO. 1 


2 


578 


yoo 


2 


91 


TKA 
{ 04 


I 1 TQ 

II /y 





80 


11 


97 


4 


5 


1 


12 


y.u4 


12 


64 


33 


2 


2 


17 


2 


32 




OiViiViJ ±ZoOio.OO-|-DZ±OOU.U 


1 


865 


ooy 


2 


20 




1 KQK 
lOoO 





66 


11 


70 


2 


4 


1 


10 


Q OK 

y.zo 


12 


28 


28 


5 


2 


33 


1 


83 


or),raa 


OiViiViJ iZoOzi .Z / +D^i / Uo.4 


1 


988 


330 


5 


65 


266 


1797 





32 


11 


91 


1 


5 


1 


30 


9.43 


12 


19 


24 


4 


1 


92 


1 


69 


or5,raa 


OiViiViJ izoDzz.DO+OziOzy . / 


2 


466 


1981 


6 


78 


1438 


1403 


1 


72 


11 


92 


10 


8 





90 


8.90 


12 


92 


45 


1 


2 


40 


2 


53 


X 


QA/fATTIO'^fiOO 1'^_Li^Oin/lc: o 
OiViiViJ iZoOzy . io-|-DZiU40.o 


1 


013 


202 





69 


176 


207 





35 


11 


70 


1 


5 


1 


53 


9.07 


12 


01 


24 


4 


2 


13 


2 


45 


V 

A. 


CTVTA/fTioQCQO CI 1 cononn i 
OiViiViJ izoDoz.Di -hoZUoUU. i 


1 


993 


1067 


17 


30 


973 


1107 


4 


52 


11 


33 


22 


7 





45 


8.45 


12 


75 


71 


5 


1 


70 


2 


46 


X,spec,mIR,,rad 


OiViiViJ iZv>Dv>4.0i -|-DziZ4i .U 


1 


219 


309 


4 


49 


289 


901 





64 


11 


66 


2 


4 


1 


06 


8.93 


12 


23 


28 


5 


2 


16 


2 


02 


SB, rad 


OiViiViJ iZoDo0.oy+Dzi4Z4. i 


2 


005 


1921 


14 


81 


1740 


1087 


7 


23 


11 


38 


37 


1 





25 


8.45 


13 


01 


71 


5 


2 


31 


2 


72 


X,spec,mIR 


oMMJ IZoOoD.YO+OZllOO.l 





557 


11 





94 


21 


24 





69 


10 


47 





8 





67 


9.36 


11 


08 


15 


2 


1 


04 


2 


44 


X,spec 


OiViiViJ J.Zo ( U / .Zl-|-0Zl4HJo. 1 


2 


484 


oOi 


4 


01 


QQft 
OOO 


yuo 





38 


11 


95 


1 


3 





92 


ft QK 

o.yo 


12 


29 


33 


2 


1 


65 


2 


09 


A.,raci 


OiViiViJ izo 1 ii.yfs-|-Dziozo. 1 


1 


992 


O / i 


2 


50 


OO i 


DOO 


1 


72 


11 


29 


8 








61 


ft dA 

O.D4 


12 


29 


45 


1 


2 


30 


2 


23 




OiViiViJ IZo 1 ±Z.UO-|-DZ1Z1Z.O 


2 


914 


1 Kf\ 


4 


96 


OAT 


DU4 





20 


12 


10 








1 


32 


y.oo 


12 


16 


24 


4 


1 


91 


2 


13 


X, spec, rad 


OiViiViJ ±Zo( ±u.Ui-|-DZUaZo.O 


2 


037 


879 


567 


07 


1091 


1399 


1 


28 


11 


93 











31 


ft TA 
o. / 4 


12 


80 


45 


1 





18 


2 


41 


X,spec 


QATAT T1 0T701 Q'T-Lfi0 1 m Q 
OiViiViJ ±Zo / Z±.0 1 -|-DZ±UoO.O 





979 


1 o 


2 


44 


91 


96 





33 


11 


44 





6 





80 


9. 77 


11 


72 


16 


9 


1 


12 


2 


49 


X,spec 


QA/fATTI'^IOm 1?_l_/10/10nQ 1 
OiViiViJ ioizUi.i f -|-4z4zUo.i 


3 


405 


4375 


37 


82 


3748 


1992 


16 


26 


11 


36 


88 


3 





18 


8. 14 


13 


34 


113 


3 


3 


08 


2 


79 


spec 


Q A T A i T 1 1 Or^ Q Q O _i_ /1 0/1 1 OO 1 
OiViiViJ ioizUo.oZ-|-4z4izy . i 


1 


544 




3 


75 


551 


560 


1 


75 


11 


50 


10 


9 





74 


8.70 


12 


51 


45 


1 


1 


99 


2 


51 


spec 


OATATTIQIOIO CO 1 /10/1/100 C: 

OiViiViJ ioiziZ.Dy+4Z44ZZ.O 


2 


805 


2095 


2 


18 


1470 


2710 


2 


22 


11 


82 


14 


4 


1 


29 


8.89 


12 


93 


38 


7 


3 


44 


2 


25 


spec, rad 


OiViiViJ loizio.z ( -|-4Zoyuu.y 


2 


565 


2361 


3387 


93 


80 


1498 





04 


12 


32 











32 


9.90 


11 


67 


15 


4 





00 


1 


24 


spec,mIR 


QA^A^HQIOOO QC:_I_/1 OQQ1 /I 1 
OiViiViJ loiZZZ.oO-|-4Zooi4.1 


2 


565 


DUU 


512 


57 


OiU 


ooy 





88 


11 


76 











31 


ft TK 
o. 1 


12 


47 


33 


2 





10 


2 


47 


spec,mIR 


OiViiViJ ioiZZt).ZU+4Z4o44.t) 


1 


038 


238 


8 


94 


252 


205 


2 


26 


11 


05 


7 


2 





36 


8.42 


12 


17 


38 


7 


1 


46 


2 


60 




OiViiViJ ioiZZo. / o+4Zoy4i.4 


1 


554 


742 


5 


47 


661 


5188 


3 


94 


11 


22 


19 


1 





59 


8.62 


12 


59 


45 


1 


2 


14 


1 


62 


rad 


CTV T TV T T 1 Q 1 OOO "JO 1 A'^AAr.A O 

OiViiViJ ioizZo.oU+4Z44o4.o 


2 


931 


1572 


16 


32 


1131 


1482 


2 


43 


11 


67 


15 


3 





61 


8.30 


12 


82 


61 


3 


1 


75 


2 


40 




OA/fATTIQIOQI 0^7 1 /lO/ICOn O 

OiViiViJ ioizoi.U f -|~4Z4DUy.U 


2 


713 


1212 


1 


42 


910 


966 


5 


13 


11 


25 


31 


8 





66 


8.70 


12 


72 


45 


1 


3 


95 


2 


49 




QA/TATTI'^TO'^O '?1_l_/10'?0/10 

OiViiViJ ±oizoz.oi-|-4Zoy4y .o 


2 


320 


2457 


8 


41 


1660 


1635 


1 


60 


12 


02 


10 


8 





87 


Q AT 

0.4 r 


12 


99 


61 


3 


2 


40 


2 


52 


mllx 


QA/TA/TTI'ilOQO 1/l_l_/10/11^f^ T 
OiViiViJ ioizoy . i4-|-4z4iOO. i 


2 


242 


1271 


3 


75 


1068 


795 


3 


32 


11 


51 


17 


9 





80 


9.01 


12 


79 


38 


7 


2 


87 


2 


64 




QATAT T1 /1 1 7/1 1 Q 1 _l_ i^OO QOQ n 
OiViiViJ 141 / 41.ol-|-OZZoZo.U 


1 


150 


948 


13 


50 


1 (U 


OTA 

z r 4 





75 


12 


01 


3 


8 





72 


8.48 


12 


65 


45 


1 


1 


28 


2 


97 


tt:> 

mlR 


OiViiViJ 141 ( 4Z.U4-|-oZoUzt). ( 





661 


403 


11 


08 


267 


216 


2 


02 


11 


12 


13 


8 





40 


8.23 


12 


19 


45 


1 


1 


74 


2 


61 




OiViiViJ 141 / oU.OU+oZolUl.U 


2 


128 


343 


5 


31 


241 


808 


1 


91 


11 


10 


12 


5 





66 


8.58 


12 


15 


38 


7 


2 


29 


1 


99 




OiViiViJ 141C5UU.4U-|-OlZfSZU.O 


1 


913 


lyyy 


17 


45 


iuyo 


1 AOA 
14Z4 





39 


12 


44 


3 


3 





80 


ft fli 


12 


80 


52 


6 





89 


2 


40 




SMMJ141802. 87+523011.1 


2 


127 


1325 


10 


95 


1127 


552 


11 


82 


10 


98 


63 


5 





74 


8.34 


12 


82 


61 


3 


2 


87 


2 


83 




SMMJ141809. 00+522803. 8 


2 


712 


1213 


6 


24 


551 


1641 





33 


12 


23 


3 


3 





83 


8.64 


12 


51 


45 


1 


1 


09 


2 


04 




SMMJ141813. 54+522923. 4 


3 


484 


243 


6 


65 


240 


3967 


1 


16 


11 


31 


4 


2 





60 


8.86 


12 


15 


33 


2 


1 


35 


1 


30 


spec, rad 


SMMJ163627. 94+405811. 2 


3 


180 


1332 


18 


32 


1095 


3210 


6 


32 


11 


24 


34 








65 


8.82 


12 


80 


45 


1 


2 


41 


2 


05 


spec 


SMMJ163631. 47+405546. 9 


2 


283 


1325 


10 


28 


865 


1646 


1 


77 


11 


69 


12 


5 





96 


8.94 


12 


70 


38 


7 


2 


62 


2 


24 


spec 


SMMJ163639. 01+405635. 9 


1 


495 


250 


7 


51 


1101 


1002 


6 


49 


11 


23 


4 


8 





53 


8.94 


12 


81 


32 


7 


1 


46 


2 


56 


mIR,rad 


SMMJ163650.43+405734.5 


2 


378 


1147 


73 


43 


1191 


4030 


3 


24 


11 


57 


10 


1 





33 


8.93 


12 


84 


45 


1 


1 


12 


1 


99 


spec,mIR,rad 


SMMJ163658. 19+410523.8 


2 


454 


1769 


4 


97 


1485 


1801 


2 


55 


11 


77 


13 


8 





80 


9.04 


12 


94 


45 


1 


2 


45 


2 


43 





Table A. 3. continued. 









SFR {Mq 






SSFR 


log M, 


A/burst /A/. 


M,/Lk 


log Ma 


log LiR 


Td 


Av 






SMG 




SED 


UV 


IR 


radio 


(Gyr^^) 


(Mq) 


(%) 


(Afm/Lo) 


(Mq) 


(Lq) 


(K) 


(mag) 


q 


AGN? 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


SMMJ163658. 78+405728.1 


1.190 


129 


8.41 


171 


274 


0.56 


11.49 


0.9 


0.59 


9.13 


11.10 


24.4 


1.14 


2.31 




SMMJ163704. 34+410530. 3 


0.840 


205 


2.50 


157 


74 


2.22 


10.85 


13.1 


0.92 


9.13 


11.96 


33.2 


2.43 


2.84 




SMMJ163706. 51+405313. 8 


2.374 


2020 


9.63 


1478 


1344 


1.83 


11.91 


10.9 


0.76 


9.06 


12.94 


45.1 


1.99 


2.56 


spec 


SMMJ221724. 69+001242.1 


0.510 


154 


14.71 


43 


62 


0.21 


11.31 


3.4 


0.63 


9.36 


11.40 


15.9 


0.69 


2.36 




SMMJ221725. 97+001238. 9 


3.094 


2499 


2.29 


1935 


1353 


2.10 


11.96 


12.5 


1.58 


9.38 


13.05 


38.7 


3.27 


2.67 




SMMJ221733. 02+000906.0 


0.926 


447 


0.72 


381 


333 


0.88 


11.64 


4.5 


1.56 


9.16 


12.35 


33.2 


2.58 


2.58 




SMMJ221733. 12+001120. 2 


0.652 


31 


2.77 


59 


62 


0.28 


11.33 


4.7 


1.79 


9.13 


11.54 


21.8 


1.27 


2.49 


rad 


SMMJ221733. 91+001352.1 


2.555 


875 


17.73 


731 


954 


0.78 


11.97 


3.8 


0.79 


9.01 


12.63 


38.7 


1.69 


2.40 




SMMJ221735. 15+001537. 2 


3.098 


594 


8.51 


536 


1627 


2.28 


11.37 


10.0 


0.68 


8.94 


12.49 


38.7 


1.93 


2.03 




SMMJ221735. 84+001558. 9 


3.089 


1969 


7.43 


1668 


1450 


6.49 


11.41 


35.0 


0.30 


8.35 


12.99 


71.5 


2.98 


2.58 




SMMJ221737. 39+001025.1 


2.614 


2991 


13.32 


2641 


2484 


7.05 


11.57 


37.1 


0.29 


8.47 


13.19 


71.5 


2.70 


2.54 




SMMJ221804. 42+002154. 4 


2.517 


1474 


15.50 


1239 


908 


6.30 


11.29 


33.5 


0.55 


8.85 


12.86 


52.6 


2.11 


2.65 


mIR 


SMMJ221806. 77+001245. 7 


3.623 


8825 


29.59 


7774 


11225 


20.81 


11.57 


109.9 


0.23 


8.25 


13.66 


113.3 


3.76 


2.36 


rad 


SMMJ030226. 17+000624. 5 


0.080 


3 


0.06 


2 


4 


0.07 


10.53 


0.4 


1.08 


8.28 


10.11 


11.4 


0.33 


2.24 


rad 


SMMJ030238. 62+001 106. 3 


0.276 


14 


0.04 


36 


44 


59.74 


8.78 


113.3 


0.36 


8.32 


11.32 


25.5 


4.22 


2.43 


spec, rad 


SMMJ030244. 82+000632. 3 


0.176 


6 


0.59 


11 


7 


0.37 


10.46 


0.0 


0.76 


8.22 


10.80 


20.1 


1.44 


2.69 




SMMJ123651. 76+621221. 3 


0.298 


6 


0.12 


14 


7 


1.02 


10.13 


0.9 


0.98 


8.91 


10.90 


13.3 


1.93 


2.79 




mean 


2.002 


1065 


68.35 


873 


1429 


3.51 


11.71 


15.6 


0.75 


9.02 


12.71 


40.1 


2.03 


2.34 




median 


2.148 


825 


6.78 


659 


1087 


1.72 


11.54 


7.2 


0.68 


8.91 


12.58 


38.7 


1.99 


2.40 




std dcv 


0.851 


1271 


395.38 


1066 


1731 


7.44 


0.55 


23.2 


0.36 


0.36 


0.60 


18.3 


0.95 


0.34 




mill 


0.080 


3 


0.04 


2 


4 


0.04 


8.78 


0.0 


0.18 


8.14 


10.11 


11.4 


0.00 


1.24 




max 


3.623 


8825 


3387.93 


7774 


11225 


59.74 


12.44 


113.3 


1.79 


9.90 


13.66 


113.3 


4.38 


2.97 





Note. — Column (1): SMG name. Column (2): redshift ([Chapman et al.l . l2005l ). Column (3): total star formation rate (SFR) for 
0.15 — 120 Mq stars averaged over the las t 50 M yr derived from the SED model. Column (4): SFR from UV emission interpolated 
from the SED template (using ' Kennicuttl 1 19981 ) . Column (5): SFR from IR emission (Column 12) used in all analysis throughout 
the paper (using Kcnnicutt, 199^). Column (6): SFR from radio emission derived directly from the radio data (using Bell, 200^). 
Column (7): specific SFR = SFRir/M,. Column (8): stellar mass. Column (9): Ratio of the mass of gas converted to star 
during the recent starburst episode to the total stellar mass. There are values greater than 100%, because the starburst episode 
is ongoing; 0% means that non-starburst template was adopted. Column (10): stellar mass to light ratio (luminosity at rest- 
frame K was intepolated using the best SED model). Column (11): dust mass. Column (12): total 8-1000 /im infrared luminosity. 
Column (13): dust temperature. Column (14): Average extinction Ay = 2.5 log(V-band starlight unextinguished / l/-band starlight 
observed). Column (15): FIR-r adio correlation param eter (Section I5.4.1|l . Column (16): AGN flag — X: X-ray identified AGN ; 
SB: X-ray identified starburst ({Alexander et all l2005l ): spec: spectroscopically identified AGN or QSO ijChapman et al.l . l2005l ): 
mIR: mid-IR identified AGN (Section 15. 4. 2|l : rad: radio datapoint is more than 3a above the starburst model (Section 15. 4. 2|l . The 
horizontal line divides the z > 0.5 and z < 0.5 samples. This table is available in a machine-readable form in the electronic edition 
of the Journal. 
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10° 10^ 10^ 10^ 10^ 10^ 10° 10^ 10^ 10^ 10^ 10^ 
Rest Frame Wavelength (|im) 



Fig. A 1 Spectral energy distrib utions (SEDs) of SMGs. Solid lines: the best GRASIL fits. Dashed lines: SEDs of GRB 
hosts (jMichalowski et all . I2008D shown for comparison. Squares: detections with errors, in most cases, smaller than the 
size of the symbols. Arrows: Su upper limit (values mark ed at the base) . In the c ases where our fits strongly underpredict 
the observed data at 850 /im, we adopted Lir and Td of ! Chapman et all (120051) . 



Michalowski et al.: Evolution of submillimeter galaxies 



21 
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Rest Frame Wavelength (|im) 



Fig. A.l (continued). 
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Fig. A.l (continued). 
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Fig. A.l (continued). 
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Fig. A.l (continued). 
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Fig. A.l (continued). 
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Fig. A.l (continued). 
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Fig. A.l (continued). 
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Fig. A.l (continued). 
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Fig. A.l (continued). 
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Table A. 4. Compilation of star formation rate density determinations in Mq yr~^ Mpc~^ 

jOr Reference 





z 


Az 


PSFR 


-ApsFR 


+ ApsFR 


Est 


3 


780 





340 





1690 





0218 





0250 


UV 


4 


920 





330 





1089 





0300 





0414 


uv 


5 


740 





360 





1194 





0423 





0655 


UV 





350 





150 





0356 





0058 





0070 


uv 





800 





200 





0656 





0110 





0133 


uv 


1 


350 





250 





0925 





0259 





0361 


uv 


1 


500 





500 





1954 





0721 





1143 


uv 


2 


750 





750 





3076 





1135 





1799 


uv 


4 


000 





500 





1300 





0569 





1012 


uv 





150 





150 





0395 





0043 





0048 


uv 


3 


040 





250 





1603 





0174 





0196 


uv 


4 


130 





300 





1245 





0256 





0322 


uv 





700 





200 





0481 





0102 





0130 


uv 


1 


250 





250 





0652 





0154 





0201 


uv 





150 





150 





0428 





0125 





0176 


uv 





750 





250 





1019 





0297 





0420 


uv 


1 


250 





250 





1368 





0399 





0564 


uv 


1 


750 





250 





1062 





0310 





0438 


uv 





350 





150 





0289 





0043 





0051 


uv 





625 





125 





0542 





0091 





0110 


uv 





875 





125 





1050 





0307 





0433 


uv 


4 


850 





450 





0350 





0150 





0150 


uv 





700 





300 





0462 





0060 





0068 


uv 


1 


250 





250 





0668 





0195 





0276 


uv 





350 





150 





0495 





0258 





0272 


uv 





750 





250 





0733 





0109 





0089 


uv 


1 


500 





500 





0988 





0087 





0095 


uv 


2 


500 





500 





2113 





0273 





0313 


uv 


3 


500 





500 





0922 





0190 





0187 


uv 





250 





250 





0569 





0226 





0592 


uv 





750 





250 





0638 





0176 





0446 


uv 


1 


250 





250 





0901 





0218 





0665 


uv 


1 


750 





250 





0922 





0223 





0681 


uv 


2 


500 





500 





0556 





0213 





0503 


uv 


3 


500 





500 





0519 





0240 





0616 


uv 


4 


500 





500 





0653 





0374 





1145 


uv 


5 


500 





500 





0285 





0166 





0824 


uv 





315 





115 





0373 





0005 





0005 


uv 





540 





110 





0533 





0005 





0006 


uv 





765 





115 





0957 





0006 





0006 


uv 





990 





110 





1082 





0006 





0006 


uv 


3 


800 





350 





0891 





0097 





0109 


uv 


5 


000 





350 





0331 





0043 





0049 


uv 


5 


900 





300 





0224 





0038 





0045 


uv 


4 


000 





500 





0362 





0050 





0050 


uv 


4 


700 





500 





0300 





0175 





0175 


uv 


4 


900 





300 





0138 





0069 





0069 


uv 


4 


000 





500 





0300 





0025 





0025 


uv 


4 


700 





500 





0200 





0088 





0088 


uv 


4 


900 





300 





0088 





0044 





0044 


uv 


5 


850 





250 





0034 





0014 





0014 


uv 


1 


000 





500 





2080 





0990 





0990 


uv 


2 


000 





500 





3980 





1800 





1800 


uv 


3 


000 





500 





3220 





1600 





1600 


uv 


4 


000 





500 





0940 





0390 





0390 


uv 


5 


000 





500 





0410 





0160 





0160 


uv 


6 


000 





500 





1260 





0740 





0740 


uv 


2 


280 





330 





1778 





1407 





6733 


uv 


4 


000 





300 





1301 





0194 





0227 


uv 


4 


700 





300 





0745 





0316 





0550 


uv 



Giavalisco et al. (2004) : Hopkins (2004) 
Giavalisco et al. (20041 : Hopkins (2004) 
Giavalisco et al. (20041 : Hopkins (20041 
Wilson et al. (2002); Hopkins (20(34) 
Wilson ct al. (2002); Hopkins (2004) 
Wilson et al. (^02); Hopkins (2004) 
Massarotti et al. (20011 : Hopkins (20041 
Massarotti et al. (20011 : Hopkins (20041 
Massarotti et al. (2001); Hopkins (200^ 
S ulliva n et al. (2000); Hopki ns (2004) 
: Hopkins (20041 
: Hopkins (20041 
Hopkins (20041 
Hopkins (2004) 
'Hopkins (20^ 
^ Hopkins (20041 
Ifopkins (2004) 



Steidel et al 
Steidel et al 
Cowie et al. 
Cow ie et aL 
Trever et al. 
ConnoUv et al, 
Conno llv et al. 




Conno llv et al. (19971 : Hopkins (20041 
Lillv et al. (1996); Hopkins (2004j 
Lillv et al. (19961 : Hopkins (2004) 
Lillv et al. (19961 : Hopkins (2004) 
Iwata et al. (20031 ; van Brcu kelen et al. (20051 
Cowie et al. (1996l :~oniervil le et al. (20011 
Cowie et al. (19961 : Sonicrvil le et al. (2001) 
Sawicki et al. (1997)'; So merville et al . (2001) 
Sawicki et al. (1997); So merville et al. (2001) 
Sawicki et al. (W97); Somerville et al. (20Qlj 
Sawicki et al. (1997); Somerville ct al. (2001) 
Sawicki et al. (19971 : Somerville ct al. (2001) 



Pascarelle et al 
Pascarelle et al 
Pascarelle ct al 
Pascarelle ct al 
Pascarelle et al 
Pascarelle et al 
Pascarelle et al 
Pascarelle et al 
Mobasher et al. 
Mo bashcr ct al. 
Mobasher et al. 



(19,98) 
(1998) 
(1998) 
(1998) 
(19981 
(19981 
(19981 
(19981 
(2Q09) 
(2009) 
(20091 
(20091 
(2007) 



Somerville et al. (2001) 
Somerville et al. (2001) 
Somervill e et al. ^001) 
Somervill e et al. ^001) 
Somervill e et al. ^001) 
Somerville et al. ^001) 
Somerville et al. ^001) 
Somerville et al. (20011 



Mobas her et al, 
Bouwens et al. 
Bouwens et al. (2007) 
Bouwens et al. (20071 
Ouchi et al. (20041 
Ouchi et al. (20041 
Ouchi et al. (20041 
Ouchi et al. (20041 
Ouchi et al. (20041 
Ouchi et al. (2004) 
Stanwav et al. (2003) 
Thompson et al. (2006) 
Thompson et al. (20061 
Thompson et al. (20061 
Thompson ct al. (20(16) 
Thom pson ct al. (2006) 
Thompson et al. (20061 
Lv et al. (20091 
Yoshida et al. (20061 : Lv et al. 
Yoshida et al. (20061: Lv et al. 



(2009) 
(20091 
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z 


Az 


PSFR 


-ApsFR 


+ ApsFR 


Estinial 


2 


300 





400 





1500 





0223 





0262 


UV 





050 





050 





0126 





0026 





0033 


uv 


3 


200 





140 





1600 





0798 





1592 


UV 





330 





040 





0353 





0059 





0071 


uv 





545 





085 





0996 





0148 





0174 


uv 


1 


125 





205 





1283 





0240 





0295 


uv 


1 


750 





130 





1898 





0390 





0491 


uv 


2 


225 





145 





1407 





0364 





0491 


uv 


2 


750 





750 





1800 





0337 





0414 


uv 





900 





500 





0989 





0221 





0285 


[0 2 







025 





025 





0122 





0036 





0050 


[0 2 







200 





100 





0136 





0032 





0062 


[0 2 







300 





100 





0119 





0023 





0031 


[0 2 







400 





100 





0536 





0095 





0147 


[0 2 







500 





100 





0955 





0137 





0180 


[0 2 







600 





100 





0649 





0086 





0117 


[0 2 







700 





100 





0535 





0083 





0120 


[0 2 







800 





100 





0566 





0085 





0116 


[0 2 







900 





100 





0714 





0112 





0165 


[0 2 




1 


000 





100 





1146 





0208 





0320 


[0 2 




1 


100 





100 





0899 





0242 





0523 


[0 2 




1 


200 





100 





0859 





0286 





0859 


[0 2 







375 





125 





0197 





0033 





0034 


[0 2 







625 





125 





0594 





0174 





0171 


[0 2 







875 





125 





1396 





0814 





0817 


[0 2 







401 





Oil 





0240 





0080 





0080 


[0 3 







636 





010 





0720 





0160 





0160 


[0 3 







881 





014 





1070 





0350 





0350 


[0 2 




1 


193 





018 





2280 





0550 





0550 


2 




2 


750 





750 





2773 





0810 





1144 


H/3 





025 





025 





0249 





0056 





0072 


Ha 





800 





300 





1172 





0262 





0338 


Ha 


2 


200 





050 





2655 





0641 





0845 


Ha 


1 


250 





550 





2350 





0137 





0145 


Ha 





150 





150 





0151 





0020 





0022 


Ha 





900 





100 





1067 





0294 





0440 


Ha 


1 


300 





600 





2799 





0676 





0891 


Ha 





200 





100 





0324 





0042 





0048 


Ha 





022 





022 





0126 





0046 





0074 


Ha 





043 





043 





0240 





0026 





0029 


Ha 





243 





009 





0360 





0120 





0060 


Ha 





245 





007 





0240 





0060 





0060 


Ha 





100 





100 





0192 





0042 





0014 


Ha 





010 





010 





0158 





0033 





0071 


Ha 





840 





030 





1500 





0200 





0200 


Ha 


2 


230 





150 





1700 





0900 





1600 


Ha 





242 





009 





0180 





0040 





0070 


Ha 





840 





030 





1700 





0300 





0300 


Ha 


2 


300 





400 





3484 





0869 





0869 


Ha 


3 


050 





350 





2141 





0450 





0450 


Ha 





350 





150 





0365 





0169 





0314 


mid-IR 





625 





125 





0678 





0297 





0527 


mid-IR 





875 





125 





1337 





0602 





1096 


mid-IR 





215 





215 





0300 





0100 





0200 


mid-IR 





515 





085 





0700 





0100 





0200 


mid-IR 





100 





100 





0175 





0049 





0049 


mid-IR 





300 





100 





0301 





0140 





0140 


mid-IR 





300 





100 





0197 





0067 





0066 


mid-IR 





500 





100 





0349 





0088 





0088 


mid-IR 





700 





100 





0570 





0096 





0101 


mid-IR 





900 





100 





0616 





0120 





0121 


mid-IR 



Reference 



Lv et al. (2009 



Reddv et al. (2008) : 
Wvder et al. (20051 
Shim et al. (2007) 
Dahle n et al. (2007) 
Dahlcn et al. (2007) 
Dahlen et al. (2007) 
Dahlen et al. (2007) 
Dahlen et al. (2007) 
Wadadekar et al. (2006) 
Teplitz et al. (2003) : Hopkins (2004) 
Gallcg o et al. (2002) : Hopkins (2004) 
Hogg et al. (19981 : Hopkins (2004) 
Hogg et al. (1998) 
Hogg et al. (1998) 
Hogg et al. (1998) 
Hogg et al. ( 1998) 
Hogg et al. (1998) 
Hogg et al. (1998) 
Hogg et al. (1998) 
Hogg ct al. (1998) 
Hogg et al. (1998) 
Hogg et al. C1998) 
Hammer et a l 
Hammer et al 



Hopkins (,2004) 



Hopkins (2004) 
Hopkins (2004) 
Hopkins (^004) 



Hopkins (.2004) 



Hopkins C2004) 



Hopkins (2004) 
Hopkins (2004) 
Hopkins (,2004) 



Hopkins (2004) 
(1997); Hopkins (2 004) 
(1997); Hopkins ( 2004) 
Hammer et al. (1997) : Hopkins (2004) 
Hippelein et al. (2003) 
Hippelein et al. (2003) 
Hippelein et al. (2003) 
Hippelein et al. (2003) 
Pettini et al. (1998); Hopkins (2004) 
Perez-G onzalez et al. (2003); Hopkins (2004) 
Tresse et al. (;2002); Hopkins (2004) 
Moorwood et al. (2000); Hopkins (2004) 
Hopkins et al. (2000) : Hopkins (2004) 
Sullivan et al. (2000) : Hopkins (2004) 
Glazebrook et al. (1999) : Hopkins (2004) 
Yan et a l. (1999); Hopkins (2004) 
Tresse fc Maddox (1998); Hopkins (2004) 
Gallego et al. (1995) : Hopkins (2004) 
Gronwall (1999); So merville et al. (2001) 
Fuiita et al. (2003b) 
Hippelein et al. (2003) 
Brinchmann ct al. ( 2004) 
Hanish et al. (2006) 
Sobral et al. (2009) 
Geach et al. (2008) 
Shiova et al. (2008) 
ViUar et al. (20081 
Reddv et al. (2008) 
Reddv et al. (20081 
Flores et al. (19991 : Hopkins (20041 
Flores et al. (1999): Hopkins (2004) 
Flores et al. (1999) : Hopkins (2004) 
Mann et al. (2002) 
Mann et al. (20(12) 
Pozzi et al. (200^ 
Pozzi et al. (2004) 
Zheng et al. (20071 
Zheng et al. (20071 
Zheng et al. (2007) 

Zheng et al. (2007} 
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z 


Az 


PSFR 


-ApsFR 


+ ApsFR 


Estimat 


2 


300 





400 





2091 





0357 





0357 


mid-IR 


3 


050 





350 





1124 





0211 





0211 


mid-IR 


1 


000 





100 





2000 





0300 





0300 


mid-IR 


2 


000 





300 





1100 





0200 





0200 


mid-IR 





450 





150 





2750 





0220 





0220 


mid-IR 





800 





200 





4870 





0170 





0170 


mid-IR 


1 


250 





250 





7550 





0290 





0290 


mid-IR 


2 


000 





500 


1 


6590 





0580 





0580 


mid-IR 





100 





100 





0180 





0025 





0029 


mid-IR 





100 





100 





0163 





0023 





0027 


mid-IR 





100 





100 





0171 





0024 





0028 


mid-IR 





300 





100 





0384 





0019 





0020 


mid-IR 





300 





100 





0356 





0017 





0018 


mid-IR 





300 





100 





0330 





0016 





0017 


mid-IR 





500 





100 





0818 





0060 





0064 


mid-IR 





500 





100 





0927 





0089 





0098 


mid-IR 





500 





100 





0589 





0043 





0046 


mid-IR 





700 





100 





1052 





0052 





0054 


mid-IR 





700 





100 





1459 





0140 





0155 


mid-IR 





700 





100 





0951 





0047 





0049 


mid-IR 





900 





100 





1319 





0580 





1036 


mid-IR 





900 





100 





1877 





0852 





1559 


mid-IR 





900 





100 





1255 





0552 





0985 


mid-IR 


1 


200 





200 





1741 





0127 





0137 


mid-IR 


1 


200 





200 





2477 





0400 





0478 


mid-IR 


1 


200 





200 





1423 





0136 





0151 


mid-IR 


1 


600 





200 





1574 





0151 





0167 


mid-IR 


1 


600 





200 





2076 





0379 





0464 


mid-IR 


1 


600 





200 





1353 





0130 





0143 


mid-IR 


2 


000 





200 





1319 





0344 





0466 


mid-IR 


2 


000 





200 





2239 





0705 





1028 


mid-IR 


2 


000 





200 





1388 





0362 





0490 


mid-IR 


2 


400 





200 





1785 





0432 





0570 


mid-IR 


2 


400 





200 





3524 





1109 





1618 


mid-IR 


2 


400 





200 





2296 





0556 





0733 


mid-IR 





698 





618 





0102 





0014 





0014 


submm 


1 


775 





367 





0228 





0027 





0027 


submm 


2 


357 





209 





0486 





0054 





0054 


submm 


3 


101 





522 





0341 





0040 





0040 


submm 


2 


000 


1 


000 





1476 





0607 





0973 


submm 


4 


500 


1 


500 





1901 





1195 





2454 


submm 





057 





041 





0206 





0016 





0024 


submm 





138 





040 





0292 





0022 





0040 


submm 





250 





073 





0192 





0026 





0036 


submm 





454 





132 





0511 





0022 





0067 


submm 





824 





239 





0785 





0073 





0086 


submm 


2 


281 


1 


219 





1104 





0092 





0140 


submm 





005 





005 





0109 





0007 





0008 


radio 





080 





080 





0187 





0035 





0038 


radio 





010 





010 





0177 





0036 





0036 


radio 





070 





070 





0120 





0025 





0031 


radio 





206 





196 





0408 





0157 





0155 


radio 





464 





054 





0667 





0246 





0254 


radio 





623 





075 





0764 





0344 





0340 


radio 





804 





080 





1315 





0446 





0459 


radio 


1 


600 





640 





1641 





0557 





0522 


radio 





005 





005 





0209 





0000 





0000 


radio 





152 





149 





0220 





0010 





0010 


radio 





310 





210 





0331 





0074 





0076 


radio 





810 





290 





0851 





0262 





0271 


radio 


1 


500 





400 





1479 





0666 





0812 


radio 



Reference 



Reddv et al. 
Reddv et al. 



Caput: 



Caput: 



Santini et al. 



Santini et al. 



Santini et al. 



et al. 



et al. 



(2008} 
f 20081 
('20071 
('20071 
f2009l 
f2009l 

\ r2nn9i 

Santini et al. (2009) 
Perez-Go nzalez et a' 
Perez-Gon zale z et a' 
Perez-Gonzalez et a 



Perez-Gonzalez et a 



Perez-Gonzalez et a 



Perez-G onzalez et a' 
Perez-G onzalez et a. 
Perez-Gonzalez et a 



Perez-Gonzalez et a 



Perez-Gonzalez et a 



Perez-Gonzalez et a' 
jPerez-G onzalez et 
Perez-Gonzalez et a 



Perez-Gonzalez et a 



Perez-Gonzalez et a 



Perez-Gonzalez et a' 
Pere z-Gonzalez et a 
Perez-G onzalez et a. 
Perez-Gonzalez et a' 



Perez-Gonzalez et a 



Perez-G onzalez et a' 
Perez-Gonzalez et a' 



Perez-Gonzalez et a 



Perez-G onzalez et a' 
Perez-G onzalez et a' 
Perez-G onzalez et a' 
Perez-Gonzalez et a 



This work 
This work 
This work 
This work 
Barger et al. 
Barger et al. 
Pascale et al. 
Pascale et al. 
Pascale et al. 
Pascale et al. 
Pascale et a l. 
Pas cale et al. 
Condon et al, 
Sadler et al. 



(^20051 
(2005) 
(2005) 
(2005) 
(2005) 
(20051 
(2005) 
(2005) 
(20051 
(2005) 
(2005) 

(2005) 
(20051 

(2005) 

(20051 

(20^ 

(20051 

(20051 

(2005) 

(20051 

(20051 

(2005) 

(20051 

(20051 

(20051 

(2005) 

(20051 



(20001 : Hopkins (20041 
(20001 : Hopkins (20041 
(2009) 
(2009.) 
(20091 
(2009) 
(2009) 
(2009) 

^ (2002); Hopkins (2004) 
(2002); Hopkins (2004) 
Serjea nt et al. (2002); Hopkins (2004) 
Macha lski fc Go dlowski (2000); Hopkins 
Haarsma et al. (2000); Hopkins (2004) 
Haarsma et al. (2000) ; Hopkins (2004) 
Haarsma et al. (20001 ; Hopkins (20041 
Haars ma et al. (20001 : 
Haars ma et al. (20001 : 
Condon (1989); Hopkins (2004 
Mauch fc Sadler (20071 
Seymour et al. (20081 
Seymo ur et al. (20081 
Scymour_etaL (20081 



Hopkins 
Hopkins 
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Az 



PSFR 



-ApsFR +AjOsFR Estimator Reference 



Seymour ct al. f2008) 
Dumre et a l. (2009) 
Dunne et al. (2009) 
Dunne et al. (20091 
Dunne et al. (2009) 
Dunne et al. (2009) 
Dunne et al. (2009) 
Dunne et al. (2009) 
Dunne et al. (2009) 
Dunne et al. (2009) 
Dunne et al. (2009) 
Dunne et al. (2009) 
Dunne et al. (2009) 
Dunne et al. (2009) 
Dunne et al. (2009) 

Gcorgakakis ct al. (2003): Hopkins (2004) 

Madau et al. (1996) : Hopkins (200*41 
Madau et al. (1996) : Hopkins (2004) 
Madau et al. (1998); van Brcukclcn ct al. (2005) 
Madau et al. (1998); van Brcukclcn ct al. (2005) 



2, 


.450 


0, 


.550 


0, 


.1202 


0, 


.0756 


0, 


.1036 


radio 


0, 


.100 


0, 


,100 


0, 


.0087 


0, 


.0063 


0, 


,0062 


radio 


0, 


.300 


0, 


,100 


0, 


.0292 


0, 


.0107 


0, 


,0106 


radio 


0, 


.500 


0, 


,100 


0, 


.0385 


0, 


.0118 


0, 


.0106 


radio 


0, 


.700 


0. 


,100 


0, 


.0700 


0, 


.0175 


0, 


.0164 


radio 


0, 


.900 


0. 


,100 


0, 


.0781 


0, 


.0202 


0, 


.0172 


radio 


1, 


.100 


0, 


,100 


0, 


.1124 


0, 


.0262 


0, 


.0249 


radio 


1, 


.300 


0, 


,100 


0, 


.1327 


0, 


.0287 


0, 


,0257 


radio 


1, 


.500 


0, 


,100 


0, 


.2043 


0, 


.0369 


0, 


,0343 


radio 


1, 


.700 


0, 


,100 


0, 


.1788 


0, 


.0339 


0, 


,0277 


radio 


1, 


.900 


0, 


,100 


0, 


.1582 


0, 


.0256 


0, 


.0245 


radio 


2, 


.250 


0, 


,250 


0, 


.1061 


0, 


.0141 


0, 


.0125 


radio 


2, 


.750 


0, 


,250 


0, 


.1049 


0, 


.0150 


0, 


.0123 


radio 


3 


.250 


0, 


,250 


0, 


.0583 


0, 


.0078 


0, 


,0069 


radio 


4, 


.250 


0, 


,750 


0, 


.0184 


0, 


.0032 


0, 


,0026 


radio 


0, 


.150 


0, 


,150 


0, 


.0383 


0, 


.0131 


0, 


,0251 


X-ray 



2.750 
4.000 
2.750 
4.000 
3.200 
2.200 
2.960 
4.130 
3.500 
2.700 
3.850 
4.700 
6.000 
5.900 
5.900 
5.900 
6.000 
3.050 
5.000 
5.900 
5.850 
2.259 
2.379 
3.110 
3.156 
3.135 
3.140 
3.400 
3.438 
3.463 
3.690 
3.700 
4.500 
4.860 
5.690 
5.700 
5.700 
5.700 
5.700 
5.735 
6.500 
6.550 
6.578 
3.000 



0.750 
0.500 
0.750 
0.500 
0.140 
0.350 
0.260 
0.260 
0.500 
0.700 
0.450 
0.200 
0.200 
0.200 
0.200 
0.200 
0.400 
0.350 
0.500 
0.300 
0.250 
0.053 
0.023 
0.020 
0.025 
0.045 
0.040 
0.030 
0.033 
0.982 
0.060 
0.220 
0.064 
0.030 
0.090 
0.100 
0.050 
0.050 
0.050 
0.062 
0.050 
0.050 
0.002 
1.000 



>0.0607 
>0.0189 
>0.0290 
>0.0110 
>0.0033 
>0.0372 
>0.0370 
>0.0161 
>0.0442 
>0.0282 
>0.0166 
>0.0147 
>0.0360 
>0.0070 
>0.0073 
>0.0221 
>0.0050 
>0.0321 
>0.0137 
>0.0003 
>0.0034 
>0.0054 
>0.0024 
>0.0120 
>0.0130 
>0.0043 
>0.0300 
>0.0060 
>0.0100 
>0.0220 
>0.0021 
>0.0004 
>0.0100 
>0.0063 
>0.0032 
>0.0012 
>0.0018 
>0.0023 
>0.0007 
>0.0005 
>0.0036 
>0.0006 
>0.0005 
>0.0818 



UV 
UV 
UV 
UV 
UV 
UV 
UV 
UV 
UV 
UV 
UV 
UV 
UV 
UV 
UV 
UV 
UV 
UV 
UV 
UV 
UV 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
Lya 
submm 



Shim et al. (2007) 

Sawicki & Thompson (2006a.b): Lv et al. (2009) 
Sawicki & Thompson (2006a.b): Lv et al. (2009) 
Sawicki & Thompson (2006a.b); Lv et al. (2009) 
Paltani et al. (2007) : Lv ct al. (2009) 
Bouwens et al. (2003a) 
(2003a) 
(2003a) 
(2003b) 

(2006) 
(2006J 
(2004) 



Bouwens et al 
Bouwens et al 
Bouwens et al 
Bouwens et al 
Bouwens et al 
Bouwens et al 
Bu nker ct al 
Reddv et al. (2008 
Iwata et al. (2007) 
Shimasaku et al. (2005) 
Stanwav et a l. (2003) 
Nilsson et al . (2009) 
Palunas et a l. (2004) 
Gronwall et al. (2007) 
Nilsson et al. (2007) 
Ouchi ct al., (2008) 
Ku dritzki et al. (2000) 
Hu'et al. (1998) 
Cowie & Hu (1998) 
van Brcukclcn et al. (2005) 
Ouchi et al. (2008) 

Fuiita et al. (2003a) : van Breukelen et al 
Hu et al. (1998) 

Ouchi et al. (2003) : van Breukelen et al. (2005 

Ouchi et al. (2008) 

Aiiki et al. (2003) 

Malhotra fc Rhoads (2004) 

Shimasaku et al. (2006) 

Murayama et al. (2007) 

Rhoads ct al. (2003) 

Malhotra fc Rhoads (2004) 

Taniguchi et al. (2005) 

Kodaira et a l. (2003) 

Hughes et al. (1998): Hopkins (2004) 
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Note. — Lower limits indicate value not corrected for extinction. The data with double reference were taken directly form the 
compilation given in the second reference. This table is available in a machine-readable form in the electronic edition of the Journal. 



Michalowski et al.: Evolution of submillimeter galaxies 
Table A. 5. Compilation of stellar mass density determinations in logM© Mpc~^ 





z 






p* 




+Ap* 


Reference 







698 





618 


7.12 


0.08 


0.13 


This work 




1 


775 





367 


7.18 


0.11 


0.16 


This work 




2 


357 





209 


7.61 


0.08 


0.12 


This work 




3 


101 





522 


7.28 


0.07 


0.12 


This work 







100 





100 


8.48 


0.10 


0.10 


Borch et al. ( 


2006) 





300 





100 


8.34 


0.15 


0.15 


Borch et al. ( 


2006) 





500 





100 


8.32 


0.11 


0.11 


Borch et al. ( 


2006) 





700 





100 


8.33 


0.10 


0.10 


Borch et al. ( 


2006) 





900 





100 


8.17 


0.18 


0.18 


Borch et al. (2006) 





100 





100 


8.49 


0.05 


0.04 


Rudnick et al. (2003) 


1 


120 





480 


8.14 


0.10 


0.11 


Rudnick et al. (2003) 


2 


010 





400 


7.48 


0.16 


0.12 


Rudnick et al. (2003) 


2 


800 





400 


7.49 


0.14 


0.12 


Rudnick et al. (2003) 





950 





450 


8.46 


0.07 


0.07 


mckinson et aL (2003) 


1 


700 





300 


8.06 


0.13 


0.17 


Dickinson et aL (2003) 


2 


250 





250 


7.58 


0.07 


0.11 


Dickinson et al. (2003) 


2 


750 





250 


7.52 


0.14 


0.23 


Dickinson et al. (2003) 





950 





450 


8.61 


0.07 


0.07 


Dickinson et al. (2003) 


1 


700 





300 


8.22 


0.12 


0.16 


Dickinson et al. (2003) 


2 


250 





250 


8.01 


0.08 


0.09 


Dickinson et al. (2003) 


2 


750 





250 


7.89 


0.15 


0.20 


Dickinson et al. (2003) 





950 





450 


8.52 


0.08 


0.07 


Dickinson et al. (2003) 


1 


700 





300 


7.97 


0.17 


0.17 


Dickinson et al. (2003) 


2 


250 





250 


7.36 


0.08 


0.11 


Dickinson et al. (2003) 


2 


750 





250 


7.27 


0.18 


0.27 


Dickinson et al. (2003) 





375 





125 


8.65 


0.17 


0.12 


Cohen (2002 







650 





150 


8.65 


0.01 


0.08 


Cohen (2002 







925 





125 


8.62 


0.09 


0.08 


Cohen (2002 







500 





100 


8.83 


0.04 


0.04 


Drorv et al. 


2004. 2005) 





700 





100 


8.76 


0.04 


0.04 


Drorv et al. 


2004, 2005) 





900 





100 


8.60 


0.04 


0.04 


Drorv et al. 


2004, 2005) 


1 


100 





100 


8.55 


0.04 


0.04 


Drorv et al. 


2004, 2005) 





500 





250 


8.50 


0.27 


0.27 


Drorv et al. 


2004, 2005) 





500 





250 


8.51 


0.17 


0.17 


Drorv et al. 


2004, 2005) 


1 


000 





250 


8.42 


0.12 


0.12 


Drorv et al. 


2004, 2005) 


1 


000 





250 


8.29 


0.19 


0.19 


Drorv et al. 


2004, 2005) 


1 


500 





250 


8.38 


0.16 


0.16 


Drorv et al. 


2004, 2005) 


1 


500 





250 


8.03 


0.16 


0.16 


Drorv et al. 


2004, 2005) 


2 


000 





250 


8.09 


0.19 


0.19 


Drorv et al. 


2004, 2005) 


2 


000 





250 


8.04 


0.20 


0.20 


Drorv et al. 


2004, 2005) 


2 


625 





375 


8.15 


0.19 


0.19 


Drorv et al. 


2004, 2005) 


2 


625 





375 


7.78 


0.20 


0.20 


Drorv et al. 


2004, 2005) 


3 


500 





500 


7.92 


0.17 


0.17 


Drorv et al. 


2004, 2005) 


3 


500 





500 


7.68 


0.20 


0.20 


Drorv et al. 


2004, 2005) 


4 


500 





500 


7.37 


0.26 


0.26 


Drorv et al. 


2004, 2005) 


4 


500 





500 


7.43 


0.20 


0.20 


Drorv et al. 


2004, 2005) 





900 





100 


8.18 


0.13 


0.10 


Glazebrook et al. (2004) 


1 


200 





100 


7.82 


0.13 


0.01 


Glazebrook et al. (2004) 


1 


450 





150 


8.08 


0.09 


0.08 


Glazebrook et al. (2004) 


1 


800 





200 


7.69 


0.13 


0.11 


fflazebrook et al. (2004) 





500 





100 


8.32 


0.03 


0.03 


Fontana et al. (2003, 2004. 2006) 





700 





100 


8.53 


0.02 


0.02 


Fontana et al. (2003. 2004. 2006) 



1 


900 
150 






100 
150 


8.16 
8.26 


0.03 
0.02 


0.03 
0.02 


Fontana et al. (2003. 2004. 2006) 
Fontana et al. (2003. 2004. 2006) 


1 


450 





150 


7.96 


0.03 


0.03 


Fontana et al. (2003, 2004, 2006J 


1 


800 





200 


7.90 


0.04 


0.04 


Fontana et al. (2003, 2004. 2006) 


2 


500 





500 


7.60 


0.04 


0.04 


Fontana et al. (2003. 2004. 2006) 


3 


500 





500 


7.23 


0.12 


0.12 


Fontana et al. (2003. 2004. 2006) 


4 


500 





500 


7.73 


0.12 


0.12 


Fontana et al. (2003. 2004. 2006) 


5 


500 





500 


7.84 


0.12 


0.12 


Fontana et al. (2003. 2004. 2006) 



36 

Table A. 5. continued. 



Michalowski et al.: Evolution of submillimeter galaxies 





z 


t 


Az 


p* 




+Ap. 


0, 


.450 


0. 


.250 


8.51 


0, 


.04 


0, 


.24 


0, 


.850 


0. 


,150 


8.44 


0, 


.05 


0, 


.20 


1, 


.250 


0, 


.250 


8.19 


0, 


.11 


0, 


.13 


1, 


.750 


0, 


,250 


7.86 


0, 


.24 


0, 


.24 


2, 


.250 


0, 


,250 


7.65 


0, 


.24 


0, 


.24 


0, 


.500 


0. 


,250 


8.64 


0, 


.17 


0, 


.24 


1, 


.000 


0. 


,250 


8.29 


0, 


.31 


0, 


.35 


1, 


.625 


0, 


,375 


7.87 


0, 


.28 


0, 


.35 


2, 


.250 


0, 


,250 


7.92 


0, 


.42 


0, 


.26 


2, 


.850 


0, 


,350 


7.90 


0, 


.20 


0, 


.38 


0, 


.100 


0. 


,100 


8.75 


0, 


.07 


0, 


.06 


0, 


.350 


0. 


,150 


8.53 


0, 


.13 


0, 


.09 


0. 


.625 


0. 


,125 


8.56 


0, 


.05 


0, 


.06 


0, 


.875 


0, 


,125 


8.48 


0, 


.04 


0, 


.06 


0, 


.100 


0, 


,100 


8.75 


0, 


.12 


0, 


.12 


0, 


.300 


0, 


,100 


8.61 


0, 


.06 


0, 


.06 


0, 


.500 


0. 


,100 


8.57 


0, 


.04 


0, 


.04 


0, 


.700 


0. 


,100 


8.52 


0, 


.05 


0, 


.05 


0, 


.900 


0, 


,100 


8.44 


0, 


.05 


0, 


.05 


1, 


.150 


0, 


,150 


8.35 


0, 


.05 


0, 


.05 


1, 


.450 


0, 


,150 


8.18 


0, 


.07 


0, 


.07 


1, 


.800 


0, 


,200 


8.02 


0, 


.07 


0, 


.07 


2, 


.250 


0, 


,250 


7.87 


0, 


.09 


0, 


.09 


2, 


.750 


0. 


,250 


7.76 


0, 


.18 


0, 


.18 


3, 


.250 


0. 


,250 


7.63 


0, 


.14 


0, 


.14 


3, 


.750 


0, 


,250 


7.49 


0, 


.13 


0, 


.13 


0, 


.150 


0, 


,150 


8.76 


0, 


.13 


0, 


.11 


0, 


.097 


0, 


,084 


8.73 


0, 


.07 


0, 


.06 


0, 


.150 


0. 


,150 


8.72 


0, 


.01 


0, 


.01 


0, 


.950 


0. 


,450 


8.45 


0, 


.07 


0, 


.06 


1, 


.700 


0. 


,300 


7.74 


0, 


.22 


0, 


.15 


2, 


.250 


0, 


,250 


7.36 


0, 


.40 


0, 


.21 


2, 


.750 


0, 


,250 


7.30 


0, 


.92 


0, 


.27 


0, 


.500 


0, 


,250 


8.57 


0, 


.03 


0, 


.03 


1, 


.000 


0, 


,250 


8.37 


0, 


.02 


0, 


.02 


1, 


.500 


0. 


,250 


8.22 


0, 


.03 


0, 


.03 


2, 


.000 


0, 


,250 


8.10 


0, 


.04 


0, 


.04 


2, 


.500 


0, 


,250 


7.93 


0, 


.04 


0, 


.04 


3, 


.500 


0, 


,500 


7.59 


0, 


.05 


0, 


.05 


4, 


.500 


0, 


,500 


6.90 


0, 


.08 


0, 


.08 


3, 


.960 


0, 


,290 


7.01 


0, 


.06 


0, 


.05 


4, 


.790 


0. 


,250 


6.63 


0, 


.07 


0, 


.06 


6, 


.010 


0. 


,250 


6.29 


0, 


.09 


0, 


.07 


5, 


.000 


0, 


,600 


6.78 


0, 


.08 


0, 


.22 


6, 


.000 


0, 


,300 


6.40 


0, 


.00 


0, 


.51 


6, 


.000 


0, 


,500 


6.59 


0, 


.55 


0, 


.24 


0, 


.100 


0. 


,100 


8.59 


0, 


.04 


0, 


.04 


0, 


.500 


0. 


,500 


8.05 


0, 


.03 


0, 


.07 


1, 


.300 


0, 


,300 


7.87 


0, 


.04 


0, 


.07 


2, 


.000 


0, 


,400 


7.76 


0, 


.06 


0, 


.06 


2, 


.800 


0, 


,400 


7.59 


0, 


.11 


0, 


.06 


0, 


.100 


0. 


,100 


8.51 


0, 


.07 


0, 


.07 


1, 


.650 


0. 


,350 


7.91 


0, 


.15 


0, 


.02 


2, 


.500 


0. 


,500 


7.55 


0, 


.18 


0, 


.12 


3, 


.500 


0, 


,500 


7.27 


0, 


.39 


0, 


.93 


0, 


.550 


0, 


,150 


8.31 


0, 


.07 


0, 


.07 


0, 


.875 


0, 


,125 


8.30 


0, 


.10 


0, 


.10 


1, 


.200 


0. 


,200 


8.15 


0, 


.10 


0, 


.10 


0, 


.300 


0. 


,100 


8.46 


0, 


.03 


0, 


.03 


0, 


.500 


0, 


,100 


8.22 


0, 


.02 


0, 


.02 


0, 


.700 


0, 


,100 


8.25 


0, 


.02 


0, 


.02 



Reference 



Fontana c t al. f2003. 200 4. 2006) 
Fonta na et al. f2003. 2004, 200^ 
Font ana et al. l'2003. 2004. 2006) 
Font ana et al. C2003, 2004, 2006J 
Font ana et al. (2003, 2004, 2006j 
Fonta na et al. (2003, 2004, 2006j 
Fontana et al. (2003, 2004, 2006j 
Fontana et al. (2003, 2004, 20063 
Fontana et al. C2003, 2004, 2006J 
Fontana et al. (2003, 2004, 20061 
Cole et al. (2001) 
Brinchmann fc Ellis (2000) 
Brinchmann fc Ellis (2000) 



Per ez-Gonzalez ct al. 
Perez- Gonzale z et al. 
Perez-Gonzalez et al. 
Perez-Gonzalez et al. 
Perez-Gonzale z et al. 
Perez-Gonzalez et al. 



Brinchmann & Ellis (2000) 
(2008J 
(2008) 
(•20081 
(20081 
(20081 
(20081 
(20081 
(20081 
(20081 
(2008) 
(2008.) 
(20081 



Perez-Gonzalez et al. 
Perez-G onzale z et al. 
Perez- Gonzalez et al. 
Perez-Gonzalez et al. 
Perez-Gonzalez ct al. 



(20081 
' (20081 
\ f2008l 
' (20081 
(20091 



Pcrez-Gonzalcz et al 
Salucci fc Persic (1999' 
Driver et al. (20071 
Bell et al. (20031 
Conselice et al. (20051 
Conselice et al. (20051 
Cons elice et al. (20051 
Cons elice et al. (20051 
Eisne r et al. (20081 
Eisner et al. (20081 
Eisner et al. (20081 
Eisner et al 
Eisner et al 
Eisner et al 
Eisner et al 

Stark et al. 

Stark et al. (20^ 
Stark et al. (20091 
Stark et al. (20071 
Evles et al. (20071 
Yan et al. (20061 
Rudnick et al. (20(16) 
Rudnick ct al. (2006) 
Rudnick ct al. (2006) 
Rudnick et al. (20061 
Rudnick et al. (20061 
Marchesini et al . (2009) 
Marchesini et al. (2009) 
Marchesini et al. (2009) 
Marchesini e t al. (2009) 
Bundv et al. (2006 ) 
Bund v et al. (20061 
Bundv et al. (2006) 
Ilbert et al. (2010) 
Ilbert et al. (20101 
Ilbert et al. (20101 
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z 




Az 


p* 




+Ap 





900 





100 


8.32 


0.01 


0.01 


1 


100 





100 


8.09 


0.02 


0.02 


1 


350 





150 


7.93 


0.01 


0.01 


1 


750 





250 


7.72 


0.09 


0.14 





300 





100 


8.79 


0.17 


0.15 





500 





100 


8.63 


0.11 


0.12 





700 





100 


8.63 


0.10 


0.10 





900 





100 


8.73 


0.13 


0.13 


1 


100 





100 


8.54 


0.11 


0.11 


1 


350 





150 


8.43 


0.12 


0.12 


1 


750 





250 


8.17 


0.12 


0.12 





325 





225 


8.66 


0.10 


0.10 





725 





175 


8.61 


0.10 


0.10 


1 


150 





250 


8.45 


0.10 


0.10 





225 





175 


8.45 


0.01 


0.01 





550 





150 


8.34 


0.02 


0.02 





800 





100 


8.22 


0.01 


0.01 


1 


050 





150 


8.14 


0.01 


0.01 


1 


400 





200 


8.04 


0.02 


0.02 


2 


050 





450 


8.05 


0.01 


0.01 





025 





025 


8.81 


0.05 


0.05 





100 





100 


8.72 


0.03 


0.03 


1 


100 





400 


8.47 


0.11 


0.11 


1 


750 





250 


8.38 


0.21 


0.21 


2 


500 





500 


8.21 


0.14 


0.14 


4 


500 


1 


500 


7.93 


0.11 


0.11 


1 


250 





250 


8.37 


0.07 


0.06 


1 


750 





250 


8.12 


0.07 


0.06 



Reference 



Ilbert et al. (2010) 
Ilbert et a" (2010) 
Ilbert et al. (2010) 
Ilbert et al. (2010) 
Arnouts et al. (2007) 
Arnouts et al. (200r) 
Arnouts et al. (2007) 
Arnouts et al. (2007) 
Arnouts et al. (2007) 
Arnouts et al. (2007) 
Arnouts et al. (2007) 
Franceschini et al. (2006) 
Franceschini et al. (2006) 
Eranceschini et al. (2006) 
Pozzetti ct al. (2007) 
et al. (2007) 



Pozzett 



Pozzett 



Pozzett 



et al. (2007) 



et al. (2007 



Pozzetti et al. (2007) 
^zzetti ct al. {2Q^ 
Kochanek et al' (2001) 
Driver et al. (2006) 
Gwvn fc Hartwick (2005) 
Gwvn & Hartwick (2005) 
Gwvn fc Hartwick (2005) 
Gwvn fc Hartwick (2005) 
CaDuti et al. (2006) 
Caputi et al. (2006) 



Note. — This table is available in a machine-readable form in the electronic edition of the Journal. 



